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Abstract 

We extend the Kreiss-Majda theory of stabiHty of hyperbohc initial- 
boundary- value and shock problems to a class of systems, notably in- 
cluding the equations of magnetohydrodynamics (MHD), for which 
Majda's block structure condition does not hold: namely, simulta- 
neously symmetrizable systems with characteristics of variable mul- 
tiplicity, satisfying at points of variable multiplicity either a "totally 
nonglancing" or a "nonglancing and linearly splitting" condition. At 
the same time, we give a simple characterization of the block struc- 
ture condition as "geometric regularity" of characteristics, defined as 
analyticity of associated eigenprojections The totally nonglancing or 
nonglancing and linearly splitting conditions are generically satisfied 
in the simplest case of crossings of two characteristics, and likewise for 
our main physical examples of MHD or Maxwell equations for a crys- 
tal. Together with previous analyses of spectral stability carried out 
by Gardner-Kruskal and Blokhin-Trakhinin, this yields immediately a 
number of new results of nonlinear inviscid stability of shock waves in 
MHD in the cases of parallel or transverse magnetic field, and recovers 
the sole previous nonlinear result, obtained by Blokhin-Trakhinin by 
direct "dissipative integral" methods, of stability in the zero- magnetic 
field limit. Our methods apply also to the viscous case. 
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1 Introduction 



There are two cases where the analysis of hyperbohc boundary value prob- 
lems is well developed: first, when the system is symmetric and the bound- 
ary conditions are dissipative |F11 IF21 IFLj : second, when the system is hy- 
perbolic with constant multiplicities and the boundary conditions satisfy a 
Lopatinski conditions [FmTl iKrl IMa( )sl lMe3j . In both cases, the main 
energy estimate is proved using symmetrizers and integrations by parts. In 
the first case, the LP' estimate is immediate: the symmetrizer is given and 
the dissipativity property is assumed. In the second case, the construction of 
symmetrizers is delicate. It was achieved first by Kreiss (see also Chazarain- 
Piriou |ChPj ) in the strictly hyperbolic case. Kreiss' construction was next 
extended by Majda-Osher and Majda to systems satisfying a "block struc- 
ture condition", which holds in particular as soon as the eigenvalues of the 
system are semi-simple and have constant multiplicity (see |Me3j ) . Kreiss's 
symmetrizers are tangential pseudo-differential operators. For symmetriz- 
able systems with dissipative boundary conditions, a trivial, constant sym- 
metrizer is available in the form SA^, where S is the symmetrizer for the 
initial value problem, and is the normal coefficient matrix; see discussion 
above Theorem 15.61 

The constant multiplicity assumption is satisfied in many applications 
such as the linearized Euler equations or isotropic Maxwell' equation, but 
it is not satisfied in other interesting example such as MHD or Maxwell's 
equation in crystals. On the other hand, the symmetry condition is com- 
monly satisfied as a consequence of the existence of a conserved energy. 
For instance, the last two examples are symmetric. But in the symmetric 
case, the boundary conditions are not always dissipative, while the sharp 
criterion of stability is given by a Lopatinski condition (see examples in 
Majda-Osher). In particular, for shock waves considered as transmission 
problems, the Rankine-Hugoniot transmission conditions are not dissipative 
in general, and Majda's analysis of shock waves relies on the use of Kreiss 
symmetrizers. Therefore, with the goal of studying shock waves for MHD, it 
is natural to extend the construction of Kreiss symmetrizers to cases where 
the eigenvalues may have variable multiplicity. 

Such an extension is given below: for symmetric systems, we extend the 
construction of symmetrizers to a framework which contains MHD. First, 
we classify the multiple eigenvalues into categories: algebraically regular, 
geometrically regular and nonregular. We show that Majda's block structure 
condition is equivalent to geometric regularity. This indicates clearly where 
Kreiss' construction stops. Next, we extend the notion of being nonglancing 
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to multiple eigenvalues. Finally, using the symmetry of the system, we 
extend the construction of Kreiss' symetrizers to the case of eigenvalues that 
are not geometrically regular provided that they are "totally nonglancing" 
or else "linearly splitting" (defined in Section 3); see Theorem 14.21 

For applications to nonconstant coefficient or quasilinear systems, it is 
important that symmetrizers vary smoothly with respect to frequency, i.e., 
that the choice of symmetrizer be robust under changes in the various pa- 
rameters. The basic construction for the totally nonglancing case is smooth. 
However, the one in the linearly splitting case is not; thus, it is useful mainly 
in the constant coefficient case. A second possibility, under sufficiently 
strong structural assumptions, would be to perform a second microlocal- 
ization about the locus of variable multiplicity; however, we do not pursue 
this direction, lacking a motivating physical example. 

Instead, we pursue a different course, investigating in Section 6 the con- 
ditions under which one may obtain a smooth symmetrizer in the original 
frequency variables, without further microlocalization. This turns out to 
have a simple and pleasing answer, which also clarifies the relation between 
Kreiss' theory and the earlier Friedrichs theory of symmetrizable systems 
with dissipative boundary conditions. Namely, one may define at points of 
variable multiplicity an appropriate reduced system, with reduced boundary 
condition, involving only the modes that change multiplicity. Then, neces- 
sary and sufficient conditions for existence of a smooth symmetrizer are that 
the full system satisfy the uniform Lopatinski condition, and the reduced 
problem admit a constant symmetrizer, or equivalently be symmetrizable 
with maximally dissipative boundary condition; see Theorems 15 . 6M^. 101 In 
the simplest case of a double root, for which the reduced system is 2 x 2, 
it turns out that the Lopatinski conditions is equivalent to symmetrizabil- 
ity/maximal dissipativity (Appendix D; see also |MaOsj ). Thus, we find that 
existence of a smooth symmetrizer requires only the Lopatinski condition 
and nothing further; see Theorem 16. 191 

Finally, in Section 7, we apply the tools we have developed to the mo- 
tivating problem of shock stability in MHD. In this case, it turns out that 
all characteristics are at least algebraically regular. More, they are either 
geometrically regular or totally nonglancing; see calculations. Appendix A. 
Thus, our results yield existence of smooth symmetrizers, and consequent 
linearized and nonlinear stability, provided the uniform Lopatinski condi- 
tion is satisfied, thus reducing the stability problem as in the constant- 
multiplicity case to a (nontrivial!) linear algebraic calculation. This imme- 
diately gives several new results and recovers the sole existing result, ob- 
tained by Blokhin and Trakhinin |BTHlBT4] by direct methods, of stability 
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in the zero-magnetic field limit; see Corollary 17.61 

Though we concentrate here on inviscid, i.e., first-order problems, the 
methods introduced extend in straightforward fashion to the viscous case 
considered in |ZlllZ31l(;MWZllinWMZ2..GWMZ3..GMWZ4j : see Appendix 
E. In particular, this extends the results obtained in |Z2llZ3j for compressible 
Navier-Stokes equations to the case of (viscous) MHD. 

2 Multiple eigenvalues of hyperbolic systems 

In this section, we recall several basic definitions and start a classification 
of multiple eigenvalues. 

2.1 Basic definitions 

Consider an x first order system with symbol 

d 

(2.1) L(p, r, = rid + A{p, = rid + ^ ^jAj{p). 

i=i 

The characteristic polynomial is A(p, r, ^) = detL(p, r, ^). Recall the fol- 
lowing definition. 

Definition 2.1. (i) The homogenous polynomial ^{r]) is hyperbolic in the 
real direction v if and only if Tr{i') ^ and for all real r]' ^ Mz^, all the roots 
in z £ C of tt{zi' + rj') = are real. 

(a) The system (|2.1() is hyperbolic in the direction v if its characteristic 
polynomial is. 

(Hi) L is symmetric hyperbolic in the direction v in the sense of Friedrichs 
if there is a smooth self adjoint matrix S{p), called a Friedrichs symmetrizer, 
such that all the matrices S{p)Aj[p) are self adjoint and S{p)L{p,v) is def- 
inite positive. This implies hyperbolicity. 

Below, we always assume that L is hyperbolic in the direction dt = 
(1,0,..., 0). This means that all the eigenvalues of A{p, ^) are real when ^ 
is real. 

Suppose that A(p, r, ^) = 0. Then — r is an eigenvalue of A(p, ^). Its 
algebraic multiplicity m is the multiplicity of — r as a root of the polynomial 
A(p, -j^). Its geometric multiplicty is rag = dimkerL(p, r, ^). The eigen- 
value is semi-simple when rUg = m. Recall that for symmetric hyperbolic 
systems, the eigenvalues are always semi-simple. 
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The simplest cases of multiple roots are described in the following defi- 
nition. 

Definition 2.2. Consider a root {p,T,^) o/ A(p,r, ^) = 0, of algebraic 
multiplicity m in t. 

i) {p,T,(^) is algebraically regular, if on a neighborhood to of {p,^) there 
are m smooth real functions Xj{p,£,), analytic in ^, such that for {p,(,) S lo: 

m 

(2.2) ^{jp,T,i)=e{p,T,i)J{{T + \j{p,i)) 

i=i 

where e is a polynomial in r with smooth coefficients such that e{p,T, ^) ^ 0. 

a) (p,T, ^) is geometrically regular if in addition there are m smooth 
functions ej{p,^) on oj with values in C^, analytic in ^, such that 

(2.3) A{p,i)ej{p,i) = A,(p,e)e,(p,0, 
and the ei, . . . , ore linearly independent. 

Remarks 2.3. a) Simple roots (i.e. m = 1), are geometrically regular. 

b) Roots of constant multiplicity are algebraically regular. In this case 
all the \j are equal. They are geometrically regular if and only if in addition 
they are semi-simple : the m = vrig vectors Cj form a smooth basis of the 
eigenspace. 

Examples 2.4. a) For MHD, the multiple eigenvalues are algebraically 
regular, but some are not geometrically regular; see Appendix A. 

b) For Maxwell's equation in a biaxial crystal, the multiple eigenvalues 
are not algebraically regular; see Appendix B. 

Systems with only geometrically regular multiple eigenvalues play an im- 
portant role, since they provide exactly the class of boundary value problems 
satisfying the block structure condition (see Theorem 13.41 below) . To treat 
examples such as MHD or non isotropic Maxwell Equations, we have to go 
beyond this class. 

2.2 Block reduction 

Consider a root (i?, t, ^) of A(p,T, ^) = 0, of algebraic multiplicity m in r. 
When {p,Cj is close to (p, ^), A(p, •,^) has exactly m roots, counted with 
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their multiplicity, close to r and there is a smooth block reduction of A{p,^): 

(2.4) U-\p,OA{p,OU{p,0 = (f^l'^^ ^(^^^)) 

with of dimension m and A has no eigenvalue in a neighborhood of r. 
Moreover, 

(2.5) A^{p,^) = —rid if r is semisimple. 

Denoting by V [resp. W ] the m first columns of U [resp. m first rows of 
U"^], there holds, 

(2.6) Aip,OV{p,0 = V{p,OA\p,0, W{p,OAip,0 = A^ip,OWip,0 



(2.7) ^^(p, = Wip, OAip, OVip, Wip, OV{p, = Id. 

Note that when L is symmetric hyperbolic, then A^ is symmetric since we 
we can choose U{p,^) such that U~^{p,^) = U*{p,^)S{p) implying that 

(2.8) W{p,0 = V*{p,OSip), A^{p,0 = V*ip,OS{p)A{p,OVip,0. 
There holds 

(2.9) A(p, r, = e{p, r, det (rid + A\t, 0) 
with e{p,T,Q / 0. 

2.3 Examples and lineeirly splitting eigenvalues 

We consider first the example of double roots for symmetric systems. In this 
case, the reduced matrix A^ has the form 

(2.10) AV,o^Ao>.oid.(;(-« H"^]^ 

with A, a and b smooth and homogeneous of degree 1 in ^. Double roots 
occur on the set 

(2.11) M = {ip,O-aip,O = bip,^) = 0} 
and ip,^,) G M- Several cases can be considered. 
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1) a = b = 0; this is the constant multiphcity case. 

2) is a smooth manifold of codimension 1, of equation {(l){p, ^) = 0}. 
This means that a and b vanish on M, and one can factor out (p, or powers 
of in o and 6. Assume for instance that 

(2.12) a = (p^a, b = (p^b, with (a, 6) ^ (0,0) at {p,Q. 

This occurs with fc = 1 for some eigenvalues of MHD. Then the eigenvalues 
are A + (p'^Xj where the \j are the smooth and distinct eigenvalues of 




The eigenvectors of are those of A. We are in the geometrically regular 
case. 

3) M. \s & smooth manifold of codimension 2 and more precisely 

(2.13) d^a{p,Q^d^b{p,Q^Q. 

Then a and 6 can be taken as independent coordinates near (p, ^), transversal 
to M.. The eigenvalues are \:ky/a? + 5^. (j9, ^) is not an algebraically regular 
point. This situation occurs for Maxwell's equations in bi-axial crystals. 

4) We give now an example of a more degenerate situation, which occurs 
for MHD. is a manifold of codimension 2, given by the equations 

(2.14) M = {(t) = '4) = Q}, d(l)Ad'ip^O. 

Moreover, a and b vanish at the second order on A4, and there is a smooth 
function c such that + b'^ = cP. For instance, this holds when 

(2.15) a = ^^-'i/'^ 6 = 2^^, c = ^^ + 'i/'^. 

In this case the eigenvalues are smooth, equal to A it c. (p, C) is an al- 
gebraically regular point. But it is not geometrically regular, since the 
eigenvectors are 




which have no limit as (<?!>, V') (0)0). 
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These examples can be generalized to higher order roots. Suppose that 
we are given a smooth conic manifold M of codimension v and a smooth 
function A on such that for all {p,S,) £ -M, is a semi-simple 

eigenvalue of A{p,S^) of multiplicity m. Suppose that (p,^) G M and r + 
A(p,O = 0. 

1) The trivial case 1^ = corresponds to constant multiplicity. 

2a) Suppose that there are no parameter p and d = 2. By homogene- 
ity, ly is at most one. When = 1, on the unit sphere |.^| = 1, Ai consists 
of isolated points, and A^{(^) is an analytic family depending on one param- 
eter. If L is symmetric, then so does , and therefore A^ is analytically 
diagonalizable: (r, ^) is a geometrically regular characteristic point. 

2b) Suppose that v = 1. Then is given by an equation = 0, and, 
extending A outside A^, A'' — Aid vanishes when </> = 0. Suppose that 

(2.16) A^(p,0 = A(p,Old + ./.'-i(p,0 

where A{p^ ^) has distinct eigenvalues. This framework extends (|2.12j) : again 
(Zi is a geometrically regular characteristic point. 

3) Suppose that v >1 and M. is given by equations (jii = . . . = (ji^ = ^ 
with dg^i, . . . , d^(t>i, linearly independent. Then 

(2.17) #(p,0 = A(p,Old + J]0,i,(p,O- 

Extending 1)2. 13(1 . consider the case where 

(2.18) A{e) ■.= ^ejAj{p,Q is strictly hyperbolic, 

3 

meaning that for all 9 £ M'^\{0}, A{6) has only real and simple eigenvalues. 
Then, for 16*1 = 1 and {p,£,) close to {p,£,), '^9jAj{p,S^) has distinct simple 
eigenvalues Xj{p,^,0), with eigenvectors ej{p,^,0). The eigenvalues of A^ 
are 

(2.19) A,(p,e) = x{p,o + mj{p,c, 

The eigenvalues are thus simple away from A4 and, in general, (p, r, ^) is 
not algebraically regular. 

In the next section, we give an intrinsic formulation of (|2.18|) . We will 
refer to this example as the linearly splitting case. 
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2.4 The tangent systems at semi-simple multiple eigenvalues 

First, we recall the following result about multiple roots of hyperbolic poly- 
nomials. 

Proposition 2.5. Assume that 7r(T, is hyperbolic in the direction dt = 
(1, 0, . . . , 0). Consider (r, C) ^ Q and assume that t is a root of multiplicity 
m > 1 of 7r{-,Q = 0. Then 

(2.20) 7r(r + r, ^ + = IL^™^ (r, + 0(|t, ^r+') 

where 7r(™) is homogeneous of degree m in (r, ^) and hyperbolic in the direc- 
tion dt. 

This means that all terms of degree smaller that m in the Taylor expan- 
sion of vr with respect to all variables, vanish. 

Consider a system (|2.1() and a root (p, r, (^) of A, of multiplicity m. 
Denote by the associated m x m block as in (|2.4() and by W and V the 
matrices defined in 1)2. 6|) (|2.7|) . Following Proposition 12.51 we also introduce 
the m-order term in the Taylor expansion of A(p, •) at (r, that we denote 
by A. 

Proposition 2.6. Assume that — r is a semi-simple eigenvalue of A{p,^) 
of order m. Then 

(2.21) A\p,^ + = -rId + A'(0 + 0(1^1') 
with 

(2.22) A!{() = Y,^jWAj{p)V, W = Wip,e), V = V{p,0 

j 

Moreover, there is e / 0, such that 

(2.23) det (rid + ^'(^)) = eA(r, e). 

Proof. Differentiating the relations (|2.(i|) (|2.7|) and using that A{p, ^) = 
-rid, yields (tTTTl) . 

The determinant A^(p, r, ^) is a factor of A{p, r, ^), as in (jSIHI)- By 

det ((r + r)Id + #(e + O) = det (rid + ^'(o) + 0(|t, er+'). 
Comparing with (jT^ . imphes ((T^ . □ 
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We will refer to L' = dt + A'{dx) as the tangent system to L at (p,T, ^). 

Remark 2.7. The operator L' is well known in the analysis of the propa- 
gation singularities or of wave packets 

u{t, x) = e^(*^+^i)/" {ao{t, x) + eai{t, x) + . . .) 

solutions to L{p,dt,dx)u = 0. The propagation is given by the laws of 
geometric optics. The principal term oq satisfies 

(2.24) ao = Val, L'(5t, 9,)a^ = 

(see (LallJMRl ; see also [EH for (IT^ ). This means that the high 
frequency oscillations are transported by L'{dt, dx)- 

Remark 2.8. By homogeneity, A^{C) = —rid and (— t, ^) is always a char- 
acteristic direction for L'. 

Example 2.9. If r + A(^) = with A an eigenvalue of A{^) of constant 
multiplicity m near ^, then 

L' = (5t + V • dx)ld, with V = d^X{Q, 

(2.25) 1 

A(r,0 = /3(r + v.O'", [3 = —d!^ ^p^r.i). 

ml — — 

L' is the transport field at the group velocity v. In this case, (|2.24|) means 
that the amplitudes qq are transported along the rays of geometric optics, 
which are the integral curve of the transport field dt + v ■ dx- 

Example 2.10. Suppose that (r, ^) is geometrically regular. With Xj and 
Bj as in (|2.2() 1)2. 3|) . in the basis ej(^), 

L' = diag(9t + vj ■ dx), v^ = d^Xj{§), 

(2.26) ^ 1 



Example 2.11. Suppose that A has the form (|2.17|) . with d^(l)i, . . . ,u,^^u 
linearly independent. Using the notation A as in (|2.18|) . the symbol of L' is 

(2.27) L'(r,0 = (T + v-e)Id + i($e) 

with V = 9^A(p, ^) and <I> = d^4>{p,S,). <I>^ = if and only if ^ is tangent 
to Aip := {4>i{p,-) = ••• = (puiPr) = 0}. Thus, the condition ()2.18|) is 
equivalent to 

(2.28) for all ^ ^ T^Mp, A'(^) has only real and simple eigenvalues. 
This motivates the following definition. 
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Definition 2.12. A multiple root (p,T, ^) is called linearly splitting transver- 
sally to a smooth manifold M = = . . . = = 0} with . . . 
analytic in ^ and d^cpi, . . . ,d^(j)u linearly independent, if Ai contains (p,^) 
and 

i) there is a smooth real valued function A(p, analytic in ^, such 
that for all {p,(,) £ M, A(p, i^) is a semi-simple eigenvalue of A{p,^) of 
constant multiplicity m, 

ii) the condition H2.28() is satisfied. 

Remark 2.13. In this case, the manifold 

(2.29) M = {{p,r,0 ■■ ip,O^M, T = -X{p,0} 

is a smooth submanifold of the characteristic variety of L which contains 
(P)Z)0- The corresponding block 

(2.30) L\p,r,0-=rU + A\p,^) 
satisfies L}'{p,T,£) = on A4. Therefore, 

(2.31) L'(r, e) = 0, for (r, C) G Tr,^Mp_. 

This means that Lf only depends on frequency variables which are transver- 
sal to Jv[. Introduce the quotient E = R^^'^/Ti-^^A^p, and denote by w the 
projection from R^'^'^ onto E. By (|2.31|) . there is L on such that 

(2.32) L'(T,e) = L(tz7(r,0) 

Because dt is transversal to wdt / 0, and the condition H2.28() is equiv- 
alent to 

(2.33) -f'(^)) & E, is strictly hyperbolic in the direction zudt. 

This approach gives a completely intrinsic definition of linear splitting. In 
particular, it shows that one can replace dt by any hyperbolic direction. 

Alternatively, consider E such that M^^*^ = E © T^^^M-p and a basis in 
E. Identifying E and E, the operator L has the form L = ^ • Ajdyj. One 
can choose E and the basis such that dt = dyo; in this case (|2.33|) or (|2.28|) 
means that L is strictly hyperbolic in the direction dyo. For example, for 
Maxwell's equation in a bi-axial crystal, the multiplicity is two and L' is a 
2x2, strictly hyperbolic system in space dimension 2, equivalent to a 2-D 
wave equation (see Appendix B). That one space dimension (at least) is lost 
when passing from L to L' is a general fact that follows from E.emark 12.81 
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Remark 2.14. If M is of codimension one, then H2.16() holds with k = 1, 
and ()2.28|) impHes that A{p, has distinct eigenvalues. This shows that if 
(p, T, ^) is linearly splitting transversally to a smooth manifold Ai of codi- 
mension 1, then (p,r, ^) is geometrically regular. 

3 The boundary block analysis 

In this section we start the investigation of noncharacteristic boundary value 
problems for hyperbolic systems. In particular, we propose an extension of 
the definition of glancing to multiple eigenvalues with variable multiplicity. 

We consider a planar boundary. Changing notations and calling (y, x) 
the spatial components, the boundary is {x = 0}. We assume that the 
boundary is not characteristic, that is 

(3.1) det Ad ^0. 

The spatial Fourier frequency variables are denoted by ^'^d r — ^7 is 

the complex time Fourier-Laplace frequency. Given a hyperbolic system L 
(|2.1|1 . we consider 

d-l 

(3.2) G{p,0 = Ad{p)-'[{T -i^)U + Y,VjMp)) 

i=i 

with C = (r,?7,7) G 1^^^+^ We denote by M:![+^ the half space {7 > 0}. By 
homogeneity, we can restrict attention to C in the unit sphere S'^; we denote 
by S'l the open half sphere S''^n {7 > 0} and by its closure S"^n {7 > 0}. 
Denoting by A the characteristic polynomial of L, there holds 

(3.3) Aip, T - n, 77, = det Adip) det (^Id + G{p, ()) . 

The consequence of hyperbolicity, is that for 7 7^ 0, G{p, C) has no real 
eigenvalue. 

3.1 Block reduction 

The goal is to construct symmetrizers for G{p, C). This is done locally, near 
each point {p, (") with E S_^_ using a diagonal block reduction of G{p, C). 

Accordingly, consider an eigenvalue — ^ of G{p, (), ( E S'\,. Denote by m' 
its algebraic multiplicity. There are N x m' , m' x m' and m' x N matrices. 
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^f>(P;C)) C (p, C) and Wh{p,C) respectively, depending smoothly on (p, C) 
near (p, (^), and holomorphically on r — ij, such that 

(3.4) GVb = VbG\ WbG = G^Wb, 

(3.5) G^ = WbGVb, WbVb = Id. 

Moreover, the spectrum of G{p,C) is 

Hyperbolicity of L implies the following. 

Lemma 3.1. With notations as above, G^ has no real eigenvalues when 
7 > 0. 

Consider = (r, 7) 7^ with 7 = 0, and a real eigenvalue — ^ of 
G'"{p,()- Thus, {p,T,r],^) is a real root of A. Denote by m its order in r 
and by A^{p, rj, ^) the m x m matrix corresponding to this root in the block 
reduction of A{p,rj,(^). There is no explicit relation between and G^ . 
However, there holds: 

Lemma 3.2. If {p,T,rj,(^) is a real root of A, then the polynomial in ^ 
det(^ + G^{p, () has real coefficients when 7 = 0. 

Moreover, for {p,T — i7, ^) in a neighborhood of {p,T,r],(,), complex 
with respect to the second and fourth argument, there holds 

(3.6) det(eid + G^(p,C)) = e{p,C,Odet{{T-i^)ld + A^{p,v,0), 

(3.7) ker(eid + G^(p,C)) = E{p, CO^ei {{r - ij)ld + A' {p,r^,0) , 

with e(p,C,0 /O andE{p,C,0 = Wbip,C)Vip,v,0- 

If in addition —t is a semi-simple eigenvalue of A(p,rj,^), then m! >m 
and E(p, ^) is of maximal rank m. 

Proof. By hyperbolicity det A^ is real and det(^Id + G) has real coefficients 
when 7 = 0. If ^ is real, then for {p,T,rj,0) close to (p, C)) the roots close 
to ^ go by conjugated pair, implying that det(^Id + G^ {p,T,rj,0)) has real 
coefficients. 

By definition of A^ and G^ , both determinant in (|3.7p are equal to A up 
to a nonvanishing factor near {p,T,r],(^). 
Moreover, 

Vip, rj, ker ((r - i7)Id + A\p, 77, o) = ker ((r - i7)Id + A{p, ry, o) 
= ker (eid + G{p, 0) = Vb{p, ker (^d + G\p, ()) ■ 
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This implies H3.7() . If — r is semi-simple, then E := ker (rid + A{p, rj, ^)) has 
dimension m and is the image of V{p, rj, It is equal to ker(^Id + G^p, ()) 
thus contained in the invariant space which is the image of Vb{p,C)- This 
implies that m' > m and that 77) is injective on E. □ 

3.2 The block structure condition 

It was pointed out by A.Majda, |Maj| (see also |MaUsj ) that Kreiss' con- 
struction of symmetrizers can be carried out if G satisfies the following block 
structure condition. 

Definition 3.3. A matrix G{p, has the block structure property near {p, Q 
if there exists a smooth matrix V on a neighborhood of that point such that 
V~^GV = diag(Gfc) is block diagonal, with blocks G^ of size x ^k, having 
one the following properties: 

i) the spectrum of GkipX) contained in {Im// 7^ 0}. 

a) Vk = 1; Gk{p, C) is real when when 7 = 0, and d^yGkip, C) / 0; 

Hi) i/fc > 1, C) has real coefficients when 7 = 0, there is /x^ G M 

such that 



(3.8) Gfc(p,C) = M,Id + 

and the lower left hand corner of dyGk{p,C) does not vanish. 

Consider a block G^ associated to an eigenvalue — ^ of G{p,(). If the 
block is elliptic, that is if ^, then G^ satisfies property i) near {p, (^), and this 
finishes the analysis when 7 > 0. 

Theorem 3.4. Consider C = (r, r/,0) and — ^ a real eigenvalue of G{p,^). 
Then, the associated block G^{p, C) satisfies Majda's block structure condition 
on a neighborhood of{p,C) if and only if{T,rj,(,) is geometrically regular for 
L is the sense of Definition \2.^ 

The if part was proved by Kreiss for strictly hyperbolic systems. It was 
next extended to hyperbolic systems with semi-simple eigenvalues of con- 
stant multiplicity in |Me3j . This proof, which only uses the factorization 
(|2.2|) and the existence of smooth eigenvectors (|2.3|) . is easily extended to 
the geometrically regular case. Surprisingly, the condition of geometric reg- 
ularity is also necessary. We postpone the independent proof of this theorem 
to Appendix C. 



1 

■•. 

■•. ■•. 1 

... 
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3.3 Glancing modes 

We want to extend the definition of glancing to general eigenvalues. Recall 
first the case of a simple eigenvalue of A{p, ry, or of an eigenvalue of con- 
stant multiplicity m. Suppose that the real eigenvalue — ^ of G{p, () satisfies 

(3.9) T + A(p,r/,e) = 0. 

Following |Krj and |Me3j . the analysis of depends on the multiplicity of ^ 
as root of (|3.9)) . Hyperbolic modes correspond to simple roots and glancing 
modes to multiple roots. Only the later yield non trivial Jordan blocks in 
(|3.8|) . The glancing condition reads 

(3.10) d^Xip,ri,O = 0. 

Introducing the transport field (|2.25jl and the speed v = (vi,...,Vrf) = 
9^ ,f A(p, r/, ^), 1)3. 1U() is equivalent to 

(3.11) Vd = 0. 

In this case, the rays which are the integral curves of the transport fields are 
tangent to the boundary. According to Example 12. 9( the conditions (|3.1()j) or 
H3.11() are satisfied if and only if the boundary is characteristic for A^™"^ , the 
principal term of A at {p,T,r],^). This motivates the following definition. 

Definition 3.5. A root {p,T,r],(^) of A, of multiplicity m in t is nonglanc- 
ing, if and only if the m-th order term in the Taylor expansion of A{p, ■) at 
(r, r/,^) satisfies 

(3.12) M'^\dx)^0 

where dx = (0, ... ,0, 1) is the (space-time) conormal to the boundary. 

Also recall that in the constant multiplicity analysis, the two cases v^^ > 
and Vrf < are different. In the first case, the rays launched at the 
boundary {x = 0} at time t = enter the domain {x > 0} for t > 0; the 
mode is said incoming. In the second case, the rays go off the domain, and 
the mode is outgoing. This has a natural extension to general eigenvalues. 
By Proposition 12.51 A^™"^ is hyperbolic in the direction dt. Introduce the 
component of dt in A*^™^ > 0. It is an open convex cone and A^'") is 
hyperbolic in any direction in r_|_. Its dual cone = {{t, y, x) : V(t; t], ^) G 
tr + yr] + xS^ > 0} is the forward cone of propagation (see e.g. |Hor j ) . In 
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the constant multiplicity case, £_)_ = {t + v'r] + Vci$, > 0} and f _,_ = M+(l,v). 
Thus, Vrf > if and only if to one of the following two equivalent conditions 
hold: 

(3.13) dx e r+, 

(3.14) r+\{0} C {x > 0}. 

Similarly, the outgoing nonglancing condition v^^ < is equivalent to one of 
the following conditions. 

(3.15) -dx £ r+, 

(3.16) r+\{o} C{x< 0}. 

This leads to the following extension. 

Definition 3.6. A root {p,T,rj,S^) of A, of multiplicity m in t, is said to 
be totally incoming if one of the two equivalent conditions (j3.13p l|3.14jl is 
satisfied. It is totally outgoing if the equivalent conditions 1)3. 13|) (|3.14l) hold. 
It is totally nonglancing if it is either totally incoming or outgoing. 

Example 3.7. We have already seen the example of constant multiplicity 
1)3. 9|) . Consider now the case of an algebraically regular root as in Example 
Then 

j 

see (|2.26|) . Each mode can be glancing, incoming or outgoing depending on 
the sign of vj^d- According to Definition 13.51 the mode is nonglancing when 
all the Vj^d are different from 0. It is totally incoming [resp., outgoing] if all 
the Vj^d are positive [resp., negative]. This explains the terminology totally 
incoming or totally outgoing. 

When — r is a semi-simple eigenvalue of of A{p, rj, ^) of multiplicity m, we 
have introduced the tangent system L'(r, ry, ^) at (p, r, rj, ^). With notations 
as in (|2.22j) . it reads 

d-l 

(3.17) L'(r, 7?, = Tld + Y, + M , Aj = WA, {p)V. 

i=i 

Lemma 3.8. Assume that — r is a semi-simple eigenvalue of A{p,r],(^). 
Then 
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i) p,T,r],^ is nonglancing if and only if the boundary is noncharacteristic 
for L' , which means that is invertihle, 

ii) p,T,rj,£^ is totally incoming [resp., outgoing] if and only L' , is hyper- 
bolic in the direction dx normal to the boundary and dx [resp. —dx ] is in 
the same component as dt; this holds if and only if the spectrum of is 
positive [resp., negative]. 

Proof. By (|2.23() . A is the characteristic polynomial of L' , up to a nonvan- 
ishing factor. This implies i). 

Since A is hyperbolic in the time direction, the eigenvalues of are 
real. In addition, the spectrum of A^ is positive if and only if the matrix 
ald+(l — a)-A^ is invertible for all a G [0, 1], thus, by definition and convexity 
of r_)_, if and only if dx belongs to this cone. In the totally outgoing case, 
the proof is similar. □ 

When — r is semi-simple and nonglancing, there is a simple link between 
the first order Taylor expansion of the reduced symbols Tld-\- A^ {p, r/, ^) and 
G\p, C) + ^Id defined at fH]! and respectively. Similarly to (fTT[|) . we 
introduce G'(C) the first order variation of at the base point: 

(3.18) G\p,C + C) = G\p,0 +G'(C) + 0(|CP). 

Proposition 3.9. Suppose that —r is a semi-simple eigenvalue of A{p,ri,S^) 
of multiplicity m and {p,T,ri,(^) is nonglancing. Then — ^ is a semi-simple 
eigenvalue ofG{p,C) of multiplicity m . Moreover, one can choose bases such 
that 

(3.19) G'(C) + eid = {A,r' ((r - i7)Id + ^'(t?, o) • 

Proof, a) We fixp = p and we omit the parameter p in the notations below. 
There holds 

(3.20) A(r + T,r/ + r?,e + 6 =e(T,r/,OA'(T,r/,0 
with e(r, rj,^) and 

m 

A'(r,r?,e) = + ^ a„„,(7?, e)r^ am^,{ri,0 = 0{\v,^r-') • 
i=i 

Thus, 

m 

(3.21) A(r,r/,0 =e(0,0,0)(r- + J]at.7)(r/,0r^) 
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where ci~^_p {-q, r) denotes the homogeneous term of degree m — j in the 
Taylor expansion of am^j. The nonglancing condition imphes that 

(3.22) a^^H0,O=(lmr with 
Thus, 

A(r, 7?, e + = eo(0, 0, 0)atk''' + 0(r+') 

showing that ^ is a root of exact order m of A(t, •) = 0. Thus, — ^ is an 
eigenvalue of G(C) of algebraic multiplicity m. 

Since — r is semi-simple, the dimension of ker(rld + A{i], ^) is equal to 
m. Since rid + A{r],^) = A^^(^Id + G{C), this space is equal to the kernel 
of (^Id + G(C) showing that the geometric multiplicity of — ^ is also equal to 
m. Thus, — ^ is semi-simple. 

b) Introduce the splittings 

(3.23) C^ = E(C)eF(C), C^ = Eb(C)eFb(C) 

where E and F [resp. Efe and ¥b\ denote the invariant spaces of A [resp. G ] 
associated to eigenvalues near and away from — r [resp. — We have shown 
that E(r/, ^) = Efe(^). Thus, one can choose bases such that the matrices V 
and Vfe occurring in (|2.6|) and ()3.4|) respectively, satisfy 

(3.24) V = V{ri,0 = Vb{Q. 

The matrix W [resp. Wb ] vanishes on F [resp. Wb] and is equal to the 
inverse of V [resp. Vf, on E [resp. E;,]. Since — r and — ^ are semi-simple, 

Ffe(C) = Range (G(C) + eid) 
(3 25) — — ~ 

= Af Range (rid + Air], o) = A^Fir,, Q 



Since A^ = WAd(p)V is invertible by Lemma 13.81 this implies that 

(3.26) Wb := WbiC) = {Ai)~^WAd and WbV = Id. 

Indeed, the matrix in the right hand side vanishes on ¥b{C) and multiplying 
it by V_ on the right gives Id. 

Differentiating (|3.4j) (|3.5j) and using that G(C) = — ^Id, implies that 

(3.27) G'iC)=mG{pX)V. 

With ([2:22]) and (l3:2fi|l the identity (l3T9|l follows. □ 
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3.4 Linearly splitting modes 

We now study the structure of when {p,T,ri,^) is a real root of A, of 
multiphcity m > 2, Hnearly sphtting transversally to M and nonglancing. 
We use the notations of Definition 12.121 We first assume that 

(3.28) dx ^ Tr,,i:Mp. 



With (|2.22)) and (|2.28|) . this is equivalent to the condition that has dis- 
tinct eigenvalues, or, since m > 2, 

(3.29) i Mid. 

In this case, there are coordinates (?/i, . . . , %-i) in M'^^^ such that M. is 
parametrized by g = (p, r/'), with -q' = (rji, . . . , rjd-u), close to q = {p, rj'): 
is locally given by 

(3.30) rj" = cl){q), ^ = -fi{q), 

with rj" = {r]d-u+i, ■ ■ ■ ,Vd-i) and /u smooth. From now on, we work in 
such coordinates. With some abuse of notations, we write 

(3.31) X{q):=X{p,rj',cP{q),-f,{q)). 

Introduce the following notations 

/3 32) q={p,v'), fi = i' - (t){q), i = (.+n{q), 

f = T + \{q), 7 = 7, (=(-?,??, 7), 

By assumption, X{q) is a semi-simple eigenvalue of A{p,r]' ,(f){q),—ii{q)^ . 
Therefore, 

(3.33) A\p,r^,i)=^{q)ld + A{q,f,,i), with i(g,0,0)=0. 

Moreover, Proposition 13. 91 implies that ii{q) is a semi-simple eigenvalue with 
multiplicity m of G{p, \{q),r]' , 0(g'), 0)) , hence 

(3.34) G\p,C) = t^{q)M + G{qX) with (5(g,0) = 0. 

To blow up the singularity near Ai, we introduce the matrices G{q, p, C) and 
A{q, p,f,^) such that for p( and p^ small enough: 

(3.35) Giq, pO =pG{q, p, Q , i(g, pi) = pA{q, p, r?, |) . 
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Proposition 3.10. Suppose that r/, ^) is a real root of A, of multiplic- 
ity m, linearly splitting transversally to the manifold Ai, nonglancing and 
satisfying ()3.28() . Then, with the notations as above, for q in a neighborhood 
of q and pC, and small, 

(3.36) det (lid + G{q, p, C)) = e{q, pC, pi) det ((f - + A{q, p, ??, 0) , 



(3.37) ker {{fu - ij)ld + A{q, p, r/, |)) = E{q, pC, pi) ker (lid + G{q, p, ()) , 

with e{q, 0, 0) ^ 0, and E{q, 0, 0) = Id. Moreover, the polynomial in i, 
det (^Id + G{q, p, C)) , has real coefficients when 7 = 0, and A(q, p,fj,(,) has 
only simple real eigenvalues when ^ is real and (17,1^) 7^ 0. 

Proof. The identities ()3.36|) and 1)3. 37(1 follow directly from Lemma 13.21 and 
(|3.26|) . through the change of variables (|3.32|) . Similarly, that det (^Id + 
G{q, p, C)) is real when 7 = follows from Lemma 13.21 

Since the mode is linearly splitting, the condition ()2.28|) implies that 
A{q,0,fj,i) has only real and simple eigenvalues when |(57,^)| = 1 and ^ is 
real. By continuity this extends to A{q, p,fi,S,) for q close to q and p small 
and 1(57,01 = 1- Since 

A{q,p,fi,i) = -A{q,pfj,pi) = \if,,i)\A{q,p\if,,i)\,{fi,i)/\ifj,i)\), 

p 

A{q, p,f},i) has only simple real eigenvalues when ^ is real and p\{'i),Q\ is 
small enough. □ 

Next, we consider the case where 

(3.38) dx G T^,iiMp. 
Equivalently, this means that 

(3.39) = aid. 

The nonglancing hypothesis implies that a ^ 0. In particular, the mode is 
totally incoming when a > and totally outgoing when o < 0. 

As in Remark 12.131 introduce the manifold M of points {p,T,r],i) with 
{p^V^O S M. and r = A(p, r/, (^). Since Lf{f,r],£,) vanishes when {f,r],$,) G 
TT,n,^-Mp, the nonglancing hypothesis implies that A^ = L/{dx) 7^ 0, thus 

(3.40) dx ^ Tr,r,,^Mp- 
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where dx is considered in the now considered as the space-time conormal to 
the boundary. However, ()3.39|) imphes that lJ{dt + adx) = 0. By H2.18() . 
if 6 is not tangent to Mp, the dimension of k ker L'(0) is at most 1, thus 

dt + adx must be tangent to Mp. This shows that, in contrast with the 
previous case, we cannot take r and ^ as independent variables transversal 
to M.. However, (|3.4U() implies that there are coordinates (r/i, . . . in 
W^^'^ such that Ai is parametrized by {q,T), with q = {p,r]i, . . . ,rid-^~i), 
and locally given by 

(3.41) ^" = ^{q,r), ^ = -^(g,r), 

with -q" = {ijd-u+i, ■ ■ ■ , Vd-i) and </>, /i smooth, analytic in r. Denote by A^b 
the manifold {rj" = (/>((7,r)}. 

Proposition 13.91 implies that — ^ is a semi-simple eigenvalue with multi- 
plicity m of G{p,T,r],0) if {p,T,rj,^) £ M.. Thus, 

(3.42) G\p,T,rj,0)= fi{q,T)ld + G{q,T,f]) with f, = r^" - 4>{q,T) 

and G((?, r, 0) = 0. Thus, there are matrices G{q, p, r, j)) such that for pi] 
small enough: 

(3.43) G{q, T, pfj) = pG{q, p, r, f]). 

Proposition 3.11. Suppose that {p,T,r],^) is a real root of A, of multiplic- 
ity m, linearly splitting transversally to the manifold A4, nonglancing and 
satisfying (|3.28|) . Then, with the notations as above, for {q,T) in a neighbor- 
hood of{q,T) and pf} small, G(q, p,T,fi) has only real and simple eigenvalues 
when fj T^O. Moreover, d^G^{p,C) = —ia^^ld 

Proof. By 1)3. 39() . we are in a totally nonglancing framework, implying that 
^ + G'(p, r, ?7, 0) is hyperbolic in the direction dx. Taking {fj,^,) as variables 
transversal to M, (|2.14l) . implies that | + G'(p, r, ry', 77, 0) = | -\- G{q,0,T,fi) 
is strictly hyperbolic in the direction dx. This shows that G{q,0,T,fi) has 
only real and simple eigenvalues when 77 7^ 0. This is preserved for (g, r) 
close to {q,r), fj ^ and p such that \pf]\ is small enough. 
By (jCTIj) . 

d-yG\p,C) = -idrG\p,Q = -iAd^ = -ia'^ld. 

□ 
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4 The Lopatinski condition and symmetrizers 
4.1 Maximal estimates 

Consider a system L{p,dt,dy,dx) = Ad{p){dx + G{p,dt,dy)), hyperbolic in 
the direction dt such that the boundary {x = 0} is not characteristic. The 
classical plane wave analysis yields the equations 

(4.1) dxU + iG{p,C)u = f, 

with C = (''" ~ ^7)"^)) identified with (r, r/, 7). We supplement the operator 
L with boundary conditions, which read, after Laplace-Fourier transform, 

(4.2) M{p,C)uix=o = g. 

We assume that M depends smoothly on p and ( S M'^"'"^\{0} and is ho- 
mogeneous of degree in ^. We allow M to depend on ^ for two reasons: 
first, it is so in the shock problem and second, the block reduction process 
leads us to consider such cases. By homogeneity, we can restrict attention 

toCGS^ = {|CI = i;7>o}. 

Recall that the hyperbolicity condition implies that for G S'^j?., the 
eigenvalues of G{p, C) are nonreal. The solutions of {dx + iG)u = 
consist of all u{x) = e~^^^h with h G E_(p, C), where E_(p, (^) denotes the 
invariant space of G{p,C) associated to the spectrum in {Im^ < 0}. Recall 
that 

(4.3) dimE_(p,C) = iV+ 

is the number of positive eigenvalues of A^. We assume that we have the 
correct number of boundary conditions, that is 

(4.4) dimker M(p, () = = N - N+, 

the number of negative eigenvalues of A^- For simplicity, we also assume 
that M is surjective, that is that M is an x N matrix. 

Following Kreiss and Majda, we are interested in the maximal estimates 
for solutions of (|iT|) (|42|) : 

(4.5) 7ll^lli2(i,^) + |n(0)p < + 

where < means that the left hand side is smaller than C times the right 
hand side where C is independent of u, f,g and C £ 5"+. A necessary condi- 
tion is that 

Ho)\ < \M{p,c)um 
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for all u £ L'^ such that dxU + iGu = 0. Thus, a necessary condition for 
is that for all ( e there holds 

(4.6) V/iGE_(p,C), \h\<C\M{p,C)h\. 

Definition 4.1. The Lopatinski determinant associated with L, M is defined 
for gamma > as 

(4.7) Z?(p,C) = det(E_(p,C),kerM(p,C)). 

We say that L, M satisfy the uniform Lopatinski condition on to if there is 
a constant c > such that 

(4.8) y{p,C)€^c^xSi: \D{p,C)\>c 

In ()4.7|) , the determinant is obtained by taking orthonormal bases in each 
space E_ and ker M (note that dimE_+dimker M = N) and independent of 
this choice. It depends on the scalar product used in to form orthonor- 
mal bases, but the condition (|4.8|) is independent of this scalar product. 
Denoting by vr^x the orthogonal projection on kerM^, then D{p,Q is the 
determinant of the restriction of 7rjyj± to E_ (in orthonormal bases). Thus, 
using that M is an isomorphism from kerM^ to C^+, uniformly bounded 
with uniformly bounded inverse for £ Sj^, we see that the condition 

(4.9) \Dip,C)\>c, 

implies ()4.6() with C depending only on c. Conversely, if ()4.6|) holds, then 
(|4.9|) is satisfied for some c > which depends only on C. 

This shows that the uniform Lopatinski condition is necessary for 1)4. 5() . 
see |Kr) . Conversely, there holds: 



Theorem 4.2 (First Main Theorem). Suppose that L is symmetric and 
that the roots of the characteristic equations are either geometrically regular 
or totally nonglancing, or nonglancing and linearly splitting transversally to 
smooth manifolds. Then, the uniform Lopatinski condition at p implies the 
maximal estimate (|4.5|) forp in a neighborhood of p. 

This was established by Kreiss f |Krj : see also |Rall IChP] ) in the strictly 
hyperbolic case by an ingenious construction of symmetrizers. It was ex- 
tended to systems with block structure by Majda and Osher Maj} IMaOsj . 



By Theorem 13 . 41 geometric regularity implies block structure. In the totally 
nonglancing case, we extend Kreiss's construction of smooth symmetrizers 
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using the symmetry of L (Theorem I5.6|) . In the presence of nonglancing 
hnearly sphtting modes, we extend further the construction of smooth sym- 
metrizers when there is at most one hnearly sphtting mode (for each Q) corre- 
sponding to a double root (Theorem l6.19() . In the general case of nonglancing 
linearly splitting modes, we prove Theorem l4.2l using nonsmooth symmetriz- 
ers. For applications to nonconstant coefficients and nonlinear equations, 
the use of nonsmooth symmetrizers seems very delicate. In this respect. 
Theorems 15.61 and 16.191 below are more important than Theorem 14.21 

Remark 4.3. When the vector bundle E_(p, (") has a continuous extension 

to (J X st^^ for some neighborhood to of p, then D extends as a continuous 

function on wxS'^^, and the uniform Lopatinski condition holds, on a possibly 
smaller neigborhood of if and only if 

(4.10) VCG^t ^(p,C)y^O, 
or, equivalently, if and only if 

(4.11) VCG:st IE-(p,C)nkerM(p,C) = {0}. 

The continuous extendability of E_ to boundary values 7 = 0, is true for 
strictly hyperbolic systems ( |Krj 1 . or when the block structure condition 
is satisfied ( |Maj| ), or, more generally, when the nongeometrically regular 
modes are totally nonglancing (see Theorem 15. 6() . In the other cases, one 
can introduce the set of limits of sequences of vectors E_(p„, C,n) as (p„, Cn) 
tends to [p, Q) and 7„ > 0: 

(4.12) ]E_(p,C)= limsup E_(p';C')- 

(p',C')-(p,C) 

7'>0 

This is a closed cone in C^, equal to E_(p, C) if 7 > and equal to the 
continuous extension of E_ when such an extension exists. If the uniform 
estimate (|4.6|) holds for p close to p and C, € 5"!^, then 

(4.13) VC G V/i G E_(p, C), \h\ < C\M{p, C)h\ 

Conversely, this implies ()4.6|1 on a neighborhood of p, with a possibly larger 
constant C. By homogeneity, it is sufficient to check (|4.13j) for h in the unit 
sphere, and by compactness this shows that this condition is equivalent to 
the following analogue of (|4.11j) : 

(4.14) VCg^+, E_(p,C)nkerM(p,C) = {0}. 
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4.2 Symmetrizers 

We recall first the method of symmetrizers. We fix a neighborhood to of p. 

Definition 4.4. A bounded symmetrizer for G{p, C) onQ, = u;xU,Uc S"^, 
is a smooth matrix on ^l, such that, for some C, c> 0, there holds 

for all [p, 0^0., 



(4.15) S(p,C) = S*(p,C), 

(4.16) |S(p,C)|<C, 

(4.17) ImS(p,C)G(p,C) >c7ld, 

It is a Kreiss symmetrizer for (G, M) if in addition 

(4.18) S(p,C) + CM*(p,C)M(p,C) > cid. 



The symmetrizer is smooth, if it extends smoothly to uj x U C uj x . 

Taking the scalar product of the equation (|4.H) with Sn and integrating 
by parts, immediately implies the following. 

Lemma 4.5. // there exists a hounded Kreiss symmetrizer for {L, M) on 
UJ X Sf_ then the maximal estimate (|4.5|) holds for all {p, C) E c<j x S'^. 

Remark 4.6. 1)4. 18|) implies that 

(4.19) S(p,C)>cId on kerM(p,C). 
Conversely, if this inequality holds, then 

{m,h) + CiK,,M^h\^>^\h\^ 

if Ci > C + c/2 + 2CVc and C > |S| as in KWi . With C2 such that 

(4.20) IvTkerM^/ll <C2|M/l|, 

this shows that ()4.19|) imphes 1)4. 18|) with c replaced hy c/2 and C by CiC|. 
This shows that one can replace the condition (|4.18l) by (|4.19|) . 

Remark 4.7. If S{p) is a symmetrizer for L, that is S{p) = S*{p) is posi- 
tive definite and all the S{p)Aj{p) are self adjoint, then T,{p) = —S{p)A(i{p) 
is a symmetrizer for G{p, C) '■ the properties (|4.15|) (|4.16|) and 1)4. 17(1 are 
satisfied. When the property 1)4. 18(1 or equivalently (|4.19|) holds, the bound- 
ary conditions are said to be strictly dissipative. In this case, the maximal 
estimates (|4.5() are satisfied. For the theory of dissipative boundary value 
problems, we refer to [FH lF2l IFLI iRal] . 
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We recall the following general fact about symmetrizers. 

Lemma 4.8. Let T, be a symmetrizer for G. Then, for 7 > 0, (^) is 
negative definite on the stable subspace K^{p,() of G{p,C,)- 

Proof. Let v satisfy dxV + iGv = 0, f (0) = u & E„, u 7^ 0. Then, f — > 
exponentially fast as x — > +00. On the other hand, using symmetry of S, 
we have 

{d/dx){T.v,v) = 2Re{T.dxV,v) = -2Re {iSGv,v) 
''^■^^^ =2{{lmT.G)v,v) >c-f\vf >0. 

Since {v, T,v) — > and is integrable as a; ^ +00, we must have, therefore 

(4.22) (u,Sn)<-c7|b||i2(K+)<0. 

as claimed. □ 

In the analysis of |Krl Maj , there is an intermediate step in the con- 
struction of symmetrizers: one construct first a family of symmetrizers S"^, 
which satisfy 1)4. 15(1 ()4.16|) and 1)4. 17() . such that the negative cone of S'^ is 
an arbitrarily small conic neighborhood of E_. Next, one uses the uniform 
Lopatinski condition : because kerM does not intersect E_, it is contained 
in the positive cone of S'^ for f k large enough implying (|4.18|) . 

Definition 4.9. Consider a family of bounded [resp. smooth ] sym- 
metrizers for G on il.'^ = uj'^ X . It is a K- family of bounded [resp. smooth] 
symmetrizers, if there are C and projectors Yl-{p,C,) onE_(p, (") such that 

(4.23) Vk > 1, V(p,C) G |n-(p,C)l < c, 

(4.24) Vk > 1, v(p,c) G ^""^ s'^(p,c) > m(K)n;n+ -nin_. 

where lij^ = Id — n_ and tu^k) +00 as k ^ +00. 

Proposition 4.10. Suppose that T,'^ is a K-family of symmetrizers for G 
on ff^. Suppose that M satisfies the uniform Lopatinski condition. Then, 
for K large enough, S'' is a Kreiss symmetrizer on $7'' for (G, M) . 

Proof. The Lopatinski condition implies that there is a constant C such that 
(|4.6)) holds. Therefore, 

|n_/i| < C\Mh\+C\M\ \U+h\. 

Thus, there are Ci and C2 such that 

< 2|n_/i|2 + 2\u+hf < Gi\Mh\^ + G2\n+h\^ - |n_/i|^ 

Therefore, if m(K) > C2, (PTTO follows. □ 
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Remark 4.11. The choice of the projector n_ is irrelevant, as long as the 
uniform bound H4.23|) holds. If n_ is another projector on E_ satisfying 
H4.23() with a constant C, then n+n„ = and 

\u+h\ = \u+n+h\ < c\u+h\, 
|n_/i| < c(|n+/i| + |n_/i|). 

Thus, 

m{K)\U+h\^ - |n_/i|2 > (m(K)/C'2 _ 2C72)|n+/i|2 _ 2C^\U_h\^. 

Therefore, changing S** to C^^S*^ we see that 1)4. 24() for and II-i- implies 
similar estimates for E'^ and Il±, with m{K) = m[K)/C'^C'^ — 2. In particular, 
we can always choose Il-{p, C) to be the orthogonal projector on E_(p, (,)■ 

4.3 Localization and block reduction 

We collect here several remarks concerning the construction of symmetrizers. 
First, the construction is local. 

Lemma 4.12. Suppose that for all {p,C) ^ uj x S'^, there are neigborhoods 
of {p, C) o-n-d self adjoint matrices C) for {p, C) G = riTl{7 > 0} 
such that the S'' form a K -family of symmetrizers for G on ^l^. Then there 
exists neigborhoods ofuJ and K- families of symmetrizers S'' on u'^ x S'^. 
If the local T,'^ are smooth, then the global Ti^ can be chosen smooth. 

Proof. By Remark l4.111 we can assume that the local symmetrizers 
on fi^ji satisfy (|4.24|) with n_ equal to the orthogonal projector on E_. 
Relabeling the families, we can also assume that they satisfy 1)4.24(1 with 
m{K) = K. 

For all K, we can find a finite covering of uJ x by open sets ilj, such 
that a /C- family of symmetrizers exists on ^j^- Consider a a partition 

of unity 1 = Y^Xj '5'+, with Xj supported in 17^. Let 

(4.25) ^^{pX) = J2^j(0^j(p^0- 

j 

Because the covering is finite, we can choose uniform constants C and c in 
(|1?T3|) . (jnni), and and the lemma follows. □ 

There is an analogous result for Kreiss symmetrizers. 
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Lemma 4.13. Suppose that for all {p,C) G x S_^, there is a neigborhood 
0, of {p, C) and a bounded [resp. smooth] Kreiss symmetriser for (G, M) on 
17+ = 0" n {7 > 0}. Then there is a neigborhood uo' of uj and a bounded 
[resp. smooth] symmetrizer on 00' x Sj^. 

Next, we consider a smooth diagonal block reduction of G on a neigh- 
borhood O of (p, C): 

(4.26) U-\pX)G{pX)U{p,0 = block diag(G,.(z,C)). 

For instance, we can consider the distinct eigenvalues /i^, /c G {1, . . . , A:}, of 
G{p,Q, small discs Dk centered at fi^ that do not intersect each other, and 
for (p, close to {p, C), the diagonal block reduction where the spectrum of 
Gk is contained in D^. 

Equivalently, taking appropriate blocks in U and U^^, one can introduce 
smooth matrices Vk and Wk such that 

(4.27) G = Y, VkGkWk, WkVj = 5,- Jd. 

For 7 > 0, we denote by Kk-{p, () the invariant subspace of Gk associ- 
ated to eigenvalues in {Im/u < 0}. Thus, 

(4.28) E_(p, = Vk{p, C)Efc,-(p, 0- 

k 

Given symmetrizers for Gk on $7, 

(4.29) S = C/-iMiag(Sfc)C/-i = J2w^^kWk. 
is a symmetrizer for G. 

Lemma 4.14. Suppose that for all k, is a K -family of symmetrizers for 
Gk on ri", then there are K-families of symmetrizers for G on il.'^. 
If the are smooth, then E'^ can be chosen smooth. 

Proof. Consider the symmetrizers Ti'^ associated by (|4.29|) to the S^. By 
Remark 14.111 we can assume that the satisfy ()4.24|) with m{K,) = k and 
Hk,- equal to the orthogonal projection on E^,-. Therefore, for h = ^ Vkhk, 
there holds 

(s'^/i, h) = Y^i^hk, hk)>KY, in^+Zifcl' - \^k,-hk?- 

Let n_ = ^k^k-Wk- It is a projector on E_ and 

Therefore, S'' = \V-'^\-'^T,'^ satisfies ^HM with m{K) = k/\U\'^\U-'^\'^ . □ 



30 



There are no simple analogue of Lemma 14.141 fro Kreiss symmetrizers, 
since the boundary conditions M do not split, in general, into boundary 
independent conditions for each block Gk- However, we give in Section 
6 an interesting result in this direction. 

We end this section with the following elementary remark. 

Remark 4.15. If S is a symmetrizer for G, then = V^*SVfc is a sym- 
metrizer for G^, since S^G^ = V^TiGVk- In particular, if the original system 
L is hyperbolic symmetric, with symmetrizer S{p), then = —S{p)Ad{p) 
is a symmetrizer for G and C) = — C)'S'(p)^(i(p)V'fc(p, C) is a 

smooth symmetrizer for Gk{p, C)- 

5 Proof of the maximal estimates 

In this section we prove Main Theorem 14.21 Because of its special interest 
for MHD, we also point out in Theorem 15.61 below . the particular case where 
only geometrically regular and totally nonglancing modes are present. In 
this case, smooth Kreiss symmetrizers are available, while they are not in 
the more general situation of Theorem 14.21 The construction of smooth 
symmetrizers is studied in more details in the next section. 

We fix C G S^. With /i^ denoting the distinct eigenvalues of G{p,C)i 
we consider on a neighborhood Q of {p,C)j the block reduction (|4.26j) of G 
such that the spectrum of Gkip,C) is {/^jr.}- this section, give different 
constructions of K-families of symmetrizers, depending on the properties of 
the blocks G^ ■ Be denote by the dimension of the block G^ , that is the 
algebraic multiplicity of When 7 > 0, G^ipX) has no real eigenvalue, 
we denote by Kk-{pX) the invariant space of Gk associated to eigenvalues 
in {Im/i < 0} and, in accordance with (|4.3I) . we set 

(5.1) dimEfc,_(p,C) =iVfc,+ 7>0. 
5.1 Kreiss' Theorem 

We first recall Kreiss' construction with Majda's extension. 

Proposition 5.1. Suppose that Gk satisfies the block structure condition at 
{p,C)- There are spaces _|_ and E(^._, neighborhoods of (p, C) CLiT-d a 
smooth family of symmetrizers () on 0,'^ = 0,'^ f] > 0} such that: 

(5.2) C< = Efc__ e Efc_+, dimEfc,_ = Nk,+, 

(5.3) E^(p, > ^m,+I^k,+ - Rl-Hk,- 
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where n_|_ is the projector on E_|_ in the decomposition H5.2|l . 

The block structure condition covers two cases: 

1) Elliptic modes. If is not real, then Gk satisfies the block structure 
condition. ThenE^ _ = E^ [resp. E^ = {0}] if Im^Ufc < [resp. Im/ifc > 0]; 
There is a symmetrizer S such that 

(5.4) ImSfcGfe>c>0 

and S is definite negative [resp. positive] when Im/ifc < [resp. luifik > 0]; 
one chooses S'' = E [resp. S'^ = kS]. 

Elliptic modes are the only possibility when 7 7^ 0. 

2) Geometrically regular modes. If 7 = and {p,T,r], —fi^J is geo- 
metrically regular then, by Theorem I3.4( Gk satisfies the block structure 
condition. The symmetrizers constructed in |Krj (see also |ChPj l satisfy 



(5.5) lm{^^kGk)=lE^, Ek>c^ld, 

which implies 1)4. 17() . This improvement is useful for applications to variable 
coefficients: it allows one to use standard Garding's inequalities for the 
operator with symbol Ek (see piP , Mok , Me2, MeZl J. 

When the block structure condition is satisfied, the spaces Kk-{p,C) 
have limits as 7 — > (see |Kr| IChP] ). It is trivial for elliptic modes. For 
geometrically regular modes, it is proved in |MeZ2j that the existence of 
symmetrizers as in Proposition 15. II implies that 

(5.6) E_ = lim E_(p,C) 

(P,C)-(P,C) 

7>0 

From its definition, geometric regularity is a local property, and remains 
satisfied in a neighborhood of the given point. Therefore, for {p, Q in a 
neighborhood of (p, C) with 7 > 0, the K^^-ip' ,C) have limits Ej^ _(p, C) as 
{p',C') ~^ (pX) with 7' > 0. Arguing as in |MeZ2j . this implies that the 
vector bundle E^ _ has a continous extension to 7 = on a neighborhood 
of(p,C). 

Therefore, with notations as in Proposition 15.11 there is a smooth split- 
ting 

c^fe =EeEfe,_(|5,c) 
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for (p, C) close to {p,C)- With projectors associated to this decompositions, 
by continuity, there holds locally 

s^(^,C)>|n^,+nfc,+ -2n^^_nfc,_. 

This shows (^) is a smooth K-family of symmetrizers for G^. Summing 
up, we have proved: 

Corollary 5.2. Suppose that satisfies the block structure condition at 
{p,Q. Then, 

i) there is a neighborhood Q. of {p,C) such that the vector bundle 
has a continous extension to i7 H {7 > 0}; 

a) there are K-families of smooth symmetrizers C) near {p,C)- 

5.2 Totally nonglancing modes 

It remains to consider the case where 7 = and the blocks are associated 
with non-geometrically regular real roots of the characteristic equations. We 
first consider totally nonglancing roots. 

Proposition 5.3. ii) If {PtZ^VjC) totally incoming, then, for {p,C) 0, 
neighborhood of {p, (), ]E_(p, () = E(p, () 

Hi) If {p,ZjVtO totally outgoing, then, for {p,C) in a neighborhood of 
{p,C),E.{p,O~=~{0}. 

Proof. We omit the parameter p in the notation below. We have to show 
that when 7 > 0, A(t — i^, rj, ^) has m roots in ^ close to ^ in the half space 
{Im^ > 0} [resp., in {Im^ > 0}] in the totally incoming [resp., outgoing] 
case. Since A has no real roots when 7 > 0, the number of these roots is 
constant for close to C and 7 > 0, and it is sufficient to count them when 
T = r and r] = r]. 

By there holds 

m m 

(5.7) A(™) (r, 0, e) = E a^-jT^r-^ = «o H + ' 

j=0 j=l 

where ^ and the second equality is the definition of the coefficients Cj . 
The hyperbolicity of A^™) implies that all the Cj are real and the nonglancing 
condition implies that the cj do not vanish. The roots of A("^)(— ^7, 0,^) =0 
are = i'^/cj. Since 

A(r - i7, 0, ^ + li) = e(0, 0, 0)7™ ( A(-) (-i, 0, |) + 0(7)) , 
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we see that the roots of A(r — 27, 0, ^) = close to ^ are 

(5.8) e = 5 + - + o(7). 

The assumption is that all the cj have the same sign. Thus, if the cj are 
positive [resp., negative], all the roots are in {Im^ > 0} [resp., in {Im^ > 0}] 
for 7 > small. □ 

Remark 5.4. When — r is semi-simple of multiplicity m, the geometric mul- 
tiplicity of — ^ as en eigenvalue of G{p, () is m. The nonglancing hypothesis 
is that the algebraic multiplicity is also m. In this case, the proof above 
shows that the condition that all the Cj have the same sign is necessary and 
sufficient for having E_ = E or E_ = {0}. 

Proposition 5.5. Suppose that L is symmetric hyperbolic and Gk is asso- 
ciated to a totally nonglancing root. Then, there are smooth K-families of 
symmetrizers for G}^ near (p, (). 

Proof. By assumption, there is a definite positive matrix S{p) such that the 
SAj are symmetric and there is a x m matrix such that for (p, C) close 

to {p,cy. 

Ek = VkC"', VkGk = GVk 
By R.emarks 14.71 and I4.15| the symmetric matrices 

Sfc(p, C) = -Vk*{p, C)S{p)Ad{p)Vk{p, C) . 
are symmetrizers for Gk- More precisely, there holds 

(5.9) Im SfcGfe = ^V^SVk > Ck^ld, Ck > 0. 

With notations as in 1)2. 22() . denote by -A^ ^ = W f,Ad(p)V ^ the boundary 
matrix of the tangent system at the real root (t, r/, ^^) under considera- 
tion. By 1)2. 8|) . the symmetry implies that W f. = V_lS(j)), hence 

(5.10) Sfc(p,C) = -^M- 

By Lemma 13.81 A,^ ^ is definite positive [resp. negative] when the mode is 
totally incoming [resp. outgoing]. With Proposition 15.31 this implies that 

= I incoming case, 

1^ KEfc in the outgoing case. 

are K-familes of symmetrizers for Gk- Q 
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With Corollary ()5.2|1 and the results of Section 4, this implies the next 
theorem. 

Theorem 5.6 (Second Main Theorem). Suppose that L is symmetric 
hyperbolic in the sense of Friedrichs and that all the real roots {p,T,rj,^) of 
the characteristic equation are either geometrically regular or totally nonglanc- 
ing. Then, there is a neighborhood uj of p such that: 

i) there there are smooth K-families of symmetrizers S'' on uj x S'^; 

ii) the vector bundle E_(p, C) has a continuous extension to j = 0; 

Hi) if in addition the boundary conditions satisfy the uniform Lopatinski, 
then there are smooth Kreiss symmetrizers and the maximal estimates (|4.5j) 
are satisfied. 

Example 5.7. The assumptions are satisfied by the equations of MHD, 
under appropriate conditions on the parameters. Indeed, all points are either 
geometrically regular or algebraically regular and totally nonglancing; see 
Appendix A. 

Remark 5.8. In Proposition 15.51 and Theorem 15.61 the symmetry of L is 
used at only one place, to construct definite (positive or negative) sym- 
metrizers for the blocks which are not geometrically regular. Thus, the 
symmetry assumption on L can be relaxed, as soon as one can construct 
such symmetrizers. For instance, this can be done for double linearly split- 
ting eigenvalues, see Section 6. 

Remark 5.9. The symmetrizers constructed in Corollary 15.21 and Theo- 
rem EiHl are obtained by localization and block reduction. Gluing the prop- 
erties 1)5. 4() for elliptic modes and 1)5. 5|) or (|5.9|) for real roots, we see that 
the symmetrizers satisfy the following condition which implies (|4.17|) : 

(5.12) Im SG = ^ VjEjVj , 

with either Ej definite positive on the support of Vj or Ej = 'yEj with 
Ej definite positive on the support of Vj. This improvement is useful in 
applications to variables coefficients and nonlinear problems. 

5.3 Linearly splitting modes 

Next we consider a block Gk associated to a linearly splitting semi simple 
nonglancing eigenvalue: we are given a smooth manifold M passing through 
{p^Vi^jJ a-s in Definition 12.121 
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Proposition 5.10. Suppose that Gt is associated to a linearly splitting 
eigenvalue transversally to a smooth manifold, and assume that it is nonglanc- 
ing. Then, there are neighborhoods O,*^ of {p, (^) and K-families of hounded 
symmetrizers E'^ on fi^ = fi'' n {7 > 0}. 

The main difference witli tlie previous cases is tliat the symmetrizers 
T,'^ are bounded, but not necessarily smooth, that is, they may have no 
continuous extensions to 

Proof, a) We first consider the case where 

(5.13) dx ^ T^,^^Mp. 

In local coordinates (g, C) as in (|3.32|) . the block Gk has the form 

(5.14) Gk= ^i{q)ld + Gk{qX). with ^((Z, 0) = 0. 
Introduce polar coordinates : 

(5.15) C = pL p = ICI, Ce5!^, 

where is the unit sphere in M'^"*"^ where Q lives, and S'^ is the half sphere 
7 > 0. With G as in (TU?^ . there holds 

(5.16) Gk{p, C) = /"('?)Id + pGk{q, 0- 

Introduce next the matrix Ak{q, p,f},^) as in (|3.35j) . By Proposition 13.101 
and the strict hyperbolicity of A, there holds for (g, p) in a neigborhood of 
(g,0) and C E S'^: 

(5.17) det (ad + Gk{q, p, C) = e(g, pC, pO Yl{(.f-ij) + Xj{q, p, r?, |)) , 

where e{q, 0, 0) 7^ and the Xj are smooth, analytic in ^, and real when ^ 
is real. Moreover, there are smooth vectors ej{q, p,fi,£,), analytic in ^ and 
linearly independent, such that 

(5.18) {Cld+Gk{q,p,f^,fi,f))ej{q,p,fi,i) = 0, -if = -Xj{q, p,fi,0. 

By Theorem l3.4l this shows that G{q, p, () satisfies the block structure condi- 
tion near (g,0, C), for all ( £ S^. Therefore, the Kreiss construction applies 
to G; Lemmas 14.131 HTTH and Corollary 15.21 imply the following. 
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Proposition 5.11. With notations as above, there is a neighborhood a) of 
{q, 0) such that: 

i) the vector bundle of invariant spaces ¥,k-iQ,P,C) associated to eigen- 
values of Gk in {Re/i < 0}, defined for 7 > 0, has a continuous extension 
to id X S_^_; its dimension is equal to the number of positive eigenvalues of 
Ak,d! 

ii) there are smooth K-families of symmetrizers T,'^{q,p,Q for Gk on 
By (|5.16j) . the negative space of Gk{p,C) is 

(5.19) Efc,_(p,C) =]Efc,_(g,|C|, ||)- 
Since is real, it follows that 

(5.20) S^(p,C) = S^(g,|CI,^) 

is a K-family of bounded symmetrizers for Gk, for (p, C) close to {p,C) ™d 
7 >• 

b) Next we consider the case where 

(5.21) dx e Trj^^^Mp. 

By (|3.39|) . the normal matrix of the tangent system at {q,T,r],^jJ is 

(5.22) 4M = ^fcId, a^y^O. 

In local coordinates {q,T,Q as in H3.42|) . the block Gk for 7 = 0, has the 
form 

(5.23) Gk\^=o = lJ-{q,T)ld + Gk{q,T,fi), with (^^(g, r, 0) = 0. 

We switch to polar coordinates i) = pi) with f/ E S'^ and introduce Gk as in 
(|3.43() . By Proposition 13 . 1 11 the eigenvalues of G are real and simple when 
fj 0- Therefore, there exists a neigborhood to of {q,T), a constant c > 0, 
and a smooth symmetric matrix TikiQ^T, PiV) on a; x S'^, such that 

(5.24) tk{q,T,p,7))>cU, Im(SG)=0. 
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Introduce 



(5.25) 



^k{q,r,\fi\,^) when?7/0, 



Id when ry = 0. 



These matrices are uniformly bounded for (p, C) in a neighborhood of {p, C) 
and 



(5.26) 



E',(g,C)>cId, Im(S',Gfc) 1^=0 = 0. 



Symmetrizer S^(g, (^) is independent of 7 and Gk is smooth; using Proposi- 
tion this imphes 



Therefore, S^, = — aS^ is a bounded symmetrizer for Gk on a neighborhood 
of {p, C) which is positive definite [resp. negative] when a is negative [resp. 
positive]. In the first case, by 1)5.22(1 and Lemma 13.81 the mode is totally 
outgoing and by Proposition 15.31 the negative space ^k- is {0}. In the 
second case, the mode is totally incoming and Ek- = E^. Thus, = kS^ 
in the outgoing case and = S^, in the incoming case provide us with 
K-families of bounded symmetrizers for Gk ■ □ 

Remark 5.12. If the mode is totally nonglancing, the space Efc^_ is either 
([^Nk |Q|. jg t];-ivially continous at (p, C)- If the mode is not totally 
nonglancing, then the matrix A^^ ^ has positive and negative eigenvalues. In 
this case, the negative space is given by (|5.19p : its dimension is positive and 
smaller then Nk, and it is a smooth function of C/lCl- I^i general, E(p, (") is 
not continuous at (p, C) since the limit depends on the direction of C- 

The Main Theorem l4.21 is a consequence of Corollarv l5.21 ProDositions l5.5l 
and 15.111 using the general results of Section 4. 

6 Smooth Kreiss symmetrizers 

We next investigate when symmetrizers may be chosen in smooth fashion; 
this is important for the treatment of variable-coefficient or nonlinear prob- 
lems. Consider again the boundary value problem 



(5.27) 



Im(S'fcGfc) 



7 



Sfc + 0(7l(p,C)-(p,C)l)+0(7'). 



a 



(6.1) 



dxu + iG{p, C)u = /, M{p, C)u\^=o = 9, 
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for (p, C) in a neighborhood of (p, C) with 7 = and G associated to a 
hyperbohc system L as in H3.2|) . 

We suppose that we are given two invariant spaces Eo(p, C) and C) 
for G(p, C), which depend smoothly on {p, C) and such that = Eq ® Ei. 
We denote by the restriction of G to and by E^ _ the negative space 
of Gk, for 7 > 0. The negative space for G is E_ = Eq,- © Ei^_. 

Assumption 6.1. Eo,-(p, C) has a continuous extension to 7 = near 
{p, C)- Moreover, dimker M = N — dimE_ and 

(6.2) kerM(p,C)nEo,-(p,C) = {0}. 

The latter two conditions are implied by the uniform Lopatinski condi- 
tion. 

6.1 Reduced boundary value problems 

Introduce 

(6.3) Eo(p,C) = M(p,C)Eo,_(p,C). 

This is a continuous bundle on where is a neighborhood of {p, C) and 
ri_,_ = n {7 > 0}. By Assumption lO dimEo = dimEo _. Let Fi denote 
a continuous bundle such that 

(6.4) MC^ = Eo©Ei. 

Definition 6.2. The reduced boundary conditions for block Ej is 

(6.5) Mj = TTjM\^^ 

where Hj denotes the projection on Ej in the splitting ()6.4() . Likewise, we 
define reduced lopatinski determinants A-^ := det (E-^, ker iVP) and a re- 
duced uniform Lopatinski condition consisting of a uniform lower bound on 
the modulus of the reduced Lopatinski determinant. 

Remark 6.3. Note that there is a useful asymmetry in Definition 16.21 in 
that we do not require continuity of at Q; indeed, Ei^_ is mentioned 
only in the Lopatinski determinant evaluated for 7 > 0. This allows us to 
separate off discontinuous blocks. The continuous block Eg may be further 
subdivided into arbitrary many blocks on which the negative subspace is 
continuous, each with its own reduced boundary condition. 
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Remark 6.4. The choice of Fq imphes that 

(6.6) ^iM|Eo,^ = 0. 

On the other hand, ttqM is not necessarily on Ei This is another 
asymmetry in the role of Eq and Ei. It reflects that the boundary conditions 
do not decouple in the block reduction of G. 

Proposition 6.5. The uniform Lopatinski condition is satisfied on the full 
space E if and only if the reduced uniform Lopatinski conditions are satisfied 
on each Ej. 

Proof, a) If the reduced problems satisfy the uniform Lopatinski condition 
on 0+, there are constants Cj such that for all (p, Q G 

(6.7) \uj\ < Cj\MjUj\ for uj G Ej^_. 

Let u = (no, ui) £ E_. Then, using (|6.6j) . we have 

\ui\ < Ci\ttiMu\ < CilvTil \Mu\ 

\uo\ < Co|Mno| < Co\Mu\ + Co\M\ \ui\, 

implying that 

(6.8) \u\<C\Mu\ fornGE_. 

b) Conversely, if the uniform Lopatinski condition is satisfied on ri_|_, 
then, since M = ttqM on Eq,-, ()6.8() restricted to Eo,_ implies that 

|no| < C|Mno| = C|Mon° I for no G Eq 

verifying the reduced Lopatinski condition on E'^. Next, by definition of Fo, 
there is an inverse R of M from Fo to Eq,-. By continuity of Eo,-, this 
inverse is uniformly bounded on a neighborhood of {p,C)- ^or ui £ Ei._, 
let no = RttqMui S Eq,- such that n = (no,ni) G E_, ttqMu = and 
Mu = TTiMu = Mini. Then 

(6.9) Iml < C2\u\ < C2C\Mu\ = C2C\Miui\. 

where C2 is the norm of the projection from on Eq. This verifies the 
reduced Lopatinski condition on E^. □ 
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6.2 Symmetrizers 

We shall seek symmetrizers commuting with the projectors onto Ej; we call 
such a symmetrizer consistent with the decomposition E = Eq © lEi- Note 
that if E is a symmetrizer for G, then the Sj defined by restricting the her- 
mitian form {T,u,v) to Ej is a symmetrizer for Gj (see Remark I4.15() . Thus, 
E = So © Si is a symmetrizer for G, consistent with the decomposisiton. 

We have shown in Section 5 that the notion of K-families of symmetrizers 
is well adapted to a block decomposition: if there are such families S^ for 
each block, then S*^ = diag(S^) is a k-family for the full system. Such an 
analogue is not true for Kreiss symmetrizers: a K-family depends only on G, 
while the definition of a Kreiss symmetrizer also depends on the boundary 
conditions, which do no decouple, in general, in the block reduction of G. 
However, we have the following useful partial result. 

Proposition 6.6. Consider G with boundary condition M , and a specified 
decomposition Eq, Ei in a neighborhood of {p,C) satisfying Assumption \6. 11 
If there exists a local smooth Kreiss symmetrizer S respecting this decompo- 
sition then 

i) the reduced boundary value problem on satisfies the (reduced) 
uniform Lopatinski condition, 

ii) the restriction Si of S to Ei is a smooth Kreiss symmetrizer for 
the reduced boundary value problem {Gi,Mi). 

Conversely, assume ii) and 

Hi) the reduced boundary value problem on E,^ satisfies the uniform 
Lopatinski condition and there exists a smooth K-family of symmetrizers for 
the reduced problem {Go,Mq). 

Then, there exists a local smooth Kreiss symmetrizer for {G,M). 

Proof, a) Suppose that S = block-diag{So, Si} is a local smooth Kreiss 
symmetrizer for (G, M). Then there are C > and e > such that 

(6.10) {u, Su) = {uQ, Soiio) + (lii, Siui) > e\uf - C\Mu\^. 

We have recalled in Section 4 that the existence of a Kreiss symmetrizer im- 
plies that the uniform Lopatinski condition must hold, and thus, by Propo- 
sition lf).5| the uniform reduced Lopatinski condition must hold on E*^ . 

Further, for all n G Ei, there is no,~ € Eq,- such that 7roM(n^+'Uo,-) = 0. 
By Lemma 14.81 and continuity of Eq,-, (SqUo,-, ""o,-) ^ 0. Moreover, by 
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(jHiSl), M{ui + no = TTiMui = Mini. Thus, 



(ni,Sini) > (ni,Sini) + (no~,Sono-) 
= ((ni + iio _), S(ni + no -)) 
> e(|nip + |no,-P) - C|M(ni + no, ''^ 
= e|n^|2-C|MV|2, 



verifying ii). 

b) Conversely, let Sq denote a K-family of symmetrizers for (Go,Mo). 
Then, by Proposition l4.1(H for {p, C) G the uniform Lopatinski condition 
implies that 

(6.12) (no, Sgno) + Co\MoUo\^ > |noP + m{K)\uo,+ \^ 

where m{K) — > +00 as k — > +00. By assumption, there is also 5*1 satisying 

(6.13) (ni,Sini) > |ni|^ - Ci|Mini|^ 
Therefore, 

-5(no, Sono) + (ni, Smi) + Co5\Mouo\'^ + CilMmip > 

|ni|^ + 5m(K)|no,+ |^ + d\uo\'^ ■ 

Using (|6.6() , we have for n = (no , ni ) : 

2|7riMn|2 > jMiml^ - C2|no,+ P 
2|7roMn|2 > |Mono|^ - Cajnip. 

Therefore, 

5(no, Sono) + (ni, Sini) + 2Co6\iroMu\'^ + 2Ci|7riMnp > 
|ni|^ + (5|noP + (5m(K)|no,+ |^ - C3(5|nip - C2|no,+ p. 

Hence, if Csd < 1/2, 5m(K) > C2, C > 2Co-5|7ro|2 + 2Ci|7ri|2 and S = 
diag((5So,i;i): 

(n,Sn) + C|Mnp > (5|np. 

□ 

Remark 6.7. If Go satisfies the block structure condition, then i) and in) 
are equivalent. 



42 



6.3 Smooth symmetrizers for Ei 

We investigate the existence of smooth Kreiss symmetrizers for blocks Ei 
associated to nonregular eigenvalues. We have in mind the case of linearly 
splitting eigenvalues. We first give necessary and next sufficient conditions. 
We assume, as we may, that Gi has a block diagonal form 

(6.14) Clip, C) = block-diag(Gi,fc(p, C)) 

where the Gi^k are some of the invariant blocks of G{p, C), for (p, (") close to 
{p,C) associated to pairwise distinct real eigenvalues = — of G{pX)- 
We further assume that the modes {p, r, rj, ^^) are semi-simple and nonglanc- 
ing. In this case, by Proposition 13.91 the are semi-simple eigenvalues of 
G and 

(6.15) Gi,k{up,C + C) = ^,Id + G!,4C) + O(ICI'). 
Moreover, denoting by 

d-l 

l/i,k = rid + ^ rjjA^.j + (,Ak^a 
i=i 

the tangent system at {p,T,r],^^), there holds 

d-l 

(6.16) G,^k=A-i;^d{i^-il)ld + ^rijA,^j). 

i=i 

The Taylor expansion of Gi at the base point reads 

(6.17) Gi(p,C + C) = Gi(p,C) + G;(C)+0(|CP). 
Gi {p, and G'^ are block diagonal; we further introduce 

(6.18) L'l = block-diag(Li fc) 

Proposition 6.8. Suppose that S(p, C) is a smooth symmetrizer for G . Let 
Si denote the restriction ofT, to Ei. Then Si is a smooth symmetrizer for 
Gi and 

i) Si := T,i(pX) is block diagonal and satisfies Im(SiGi) = 0; 

ii) IJ^ is hyperbolic symmetric in the sense of Friedichs. 
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Proof. Since S is a symmetrizer, Im (SG) > 07, which by restriction to Ei, 
imphes that for 7 > 0: 

(6.19) ImSiGi > c7ldEi. 
Thus, by continuity, 

Im(EiGi) > 0. 

Since is block diagonal, with real entries fij, this means that for all pair 
(k, I) with k ^ I, the nondiagonal block entries ^ ^ of are such that 

is nonnegative. This implies that P = 0, hence ^ ; = when k ^ I. thus 

(6.20) Im(EiGi) = 0. 

The diagonal block Si ^ ^ is a symmetrizer for Gi^k by Remark 14.151 
Taking the first order term in C, in Im (Si^^ ^Gi^fc), implies that 

(6.21) ImSi_,^,G'i,fc(C) > C7. 

We have used that the G_i i. = Ai^Id and that the derivatives of Si ^ are 
self-adjoint. Because G'l is a linear in changing C, io —C, implies that 
ImSij, fcG'i fc(C) = when 7 = 0. Thus, the m atrices S^^^ ■= -lLi,k,kAk4 
and S_i^^A^ ,j are self-adjoint. Moreover, 1)6. 21(1 implies that S_i^^ is definite 
positive. This shows that 5i ^ is a Friedrichs symmetrizer for L'l therefore 
that 

(6.22) 5i = block-diag(Si^fc) 

is a Friedrichs symmetrizer for L'l. □ 

Proposition 6.9. Suppose that Si is a smooth symmetrizer for Gi. Then, 
it is a smooth Kreiss symmetrizer on a neighborhood of (p, (") if and only if 
the reduced boundary condition Mi(p, (") is maximal dissipative for L/i with 



symmetrizer (|6.22j) . 

Proof. If Si is a Kreiss symmetrizer for (Gi,Mi), then 
(6.23) {T.iui,ui) + C\MiUi\'^ > c|uip. 

Specializing at (p, implies that S_iA^ = —Si is maximal negative on 
kerMi(p, ^), thus that the boundary condition Mi{p,(^) is maximal dissi- 
pative for L'l with symmetrizer S_i. 

Conversely, if Mi (p, (^) is maximal dissipative for L'l with symmetrizer 
S_i, then (|6.23j) holds at {p, C), therefore on a neighborhood of that point. □ 
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Theorem 6.10 (Third Main Theorem). Suppose that L is a hyperbolic 
system and M are boundary conditions which satisfy the uniform Lopatinski 
condition. For all (p, C) o, neighborhood of{p,C), denote 6y Eo(p, C) the 
invariant space of G corresponding to eigenvalues which are either nonreal, 
or real and associated to geometrically regular or totally nonglancing modes 
of L. Denote by Ei(p, ^) the invariant space associated to the remaining 
eigenvalues of G. Then, A ssumption is satisfied and, with notations as 
in Provosition \6. M there exists a smooth Kreiss symmetrizer for {G, M) near 
{Pi C) if OLf^d only if there is a smooth symmetrizer Si for Gi such that the 
associated tangent system Lfi ()6.18|) is hyperbolic symmetric in the sense of 
Friedrichs with symmetrizer S_i H6.22() and the reduced boundary condition 
is maximally dissipative for L'^ . 

Proof. By Corollary 15.21 Propositions 15.3115.51 and Lemma 14.121 there are 
smooth K-families of symmetrizers for Go and Assumption 16.11 is satisfied. 

If S is a smooth Kreiss symmetrizer for {G, M), then by Proposition 16.61 
Si is a smooth Kreiss symmetrizer for the reduced problem (Gi,Mi) and 
Proposition 16.91 implies that the reduced boundary condition is maximally 
dissipative for L'l . 

Conversely, if Si is a smooth symmetrizer for Gi such that the reduced 
boundary condition is maximally dissipative for L'^, then Si is a Kreiss 
symmetrizer for the reduced problem (Gi,Mi) and Proposition 16.61 implies 
that there is a smooth Kreiss symmetrizer for the full system (G, M). □ 

This reduces the analysis to the construction of smooth Kreiss sym- 
metrizers for Gi. When G is associated to a symmetric hyperbolic system L, 
with symmetrizer S, then S = —SA^ is a smooth symmetrizer for G. and 
the restriction Si of S to Ei is a smooth symmetrizer on a neighborhood 
of (p, C)- Foi' this symmetrizer, the maximal dissipativity condition of the 
reduced boundary condition holds if and only if Si(p, (") is definite positive 
on ker Mi(p, C), that is if and only if there is c > such that 

(6.24) VuGE^, Mu £ MEn (SA^u.u) > c|k|^ 

where _ means evaluation at the base point {p,C)- Therefore, 

Corollary 6.11. With assumptions and notations as in Theorem \6.1(A if the 

condition (|6.24|) holds at {p,C), then there is a smooth Kreiss symmetrizer 
for {G,M) on a neighborhood of {p,C)- 

Remark 6.12. Suppose that L is symmetric and split Gi in diagonal blocks 
Gi^k as in (|6.14j) . Then, by R,emark l4.15| there are smooth symmetrizers Si^^ 
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for Gi^k, of the form Si^^ = —V^i^{SAci)Vi^k- Therefore, there are smooth 
symmetrizers for Gi of the form 

(6.25) Si = block-diag(afcSi fc), > . 

For these symmetrizers, one can state the analogues of condition H6.24() and 
Corollary 16.111 

Remark 6.13. We have seen that a necessary condition for the exis- 
tence of a smooth symmetrizer Si is that L/^ must be hyperbolic symmetric 
in the sense of Friedrichs, with symmetrizer S_i = block-diag^i , where 
S_i f. = —Si kAk]i ^1 = block-diagSi f,. Therefore, in the construction 

of a smooth Kreiss symmetrizer Si, the first step is construct a Friedrichs 
symmetrizer S_i = block-diagS^i for such that M_i is maximal dissipa- 
tive. This defines Si , and the second is to determine a smooth symmetrizer 
for Gi such that Si(p, () = Si. Since Si is block diagonal, it is sufficient to 
construct smooth symmetrizers for Gi^k such that Si_fc(p, C) = k- 

When L is symmetric hyperbolic, the L'l have Friedrichs symmetrizers 
5i ^, as in Remark 16. 121 In generic cases, these symmetrizers are unique up 
to a constant, implying that necessarily 

(6.26) 5^1 = afc^i fc, ak>0. 
Therefore, 

(6.27) Si = block-diag(afcSi J. 

In this case, if there is a smooth Kreiss symmetrizer, there is one of the form 
(|6.25jl . indicating that this form contains almost all the possible choices of 
smooth symmetrizers. 

Remark 6.14. More precisely, symmetrizers S_i are unique up to constants 
if and only if the reduced tangent systems L/^ ^ are irreducible in the sense 
that there do not exist nontrivial constant invariant subspaces of Ai kiv^O- 
For, without loss of generality taking Af ^ diagonal, so that S_i must be 
diagonal as well, and expressing S_i ^, as a block-diagonal matrix with blocks 
consisting of distinct scalar multiples of the identity, we find that symmetry 
of i^i ^-A-J ^ implies holds if and only if each A^^ ^ shares the same block- 
diagonal structure. 

In the linearly splitting case, we find by spectral separation that these 
invariant subspaces must correspond to the limit as (r/, ^) ^ of group 
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eigenprojections associated to fixed subsets of the associated eigenvalues 
Aj^(r/,^) of the nearby perturbed system. In particular, this implies that 
these group eigenprojections are continuous along rays across the origin, 
which is sufficient to give analyticity separately in coordinates rjj, ^. By 
Hartog's Theorem, therefore, they are jointly analytic in and thus 

there is an analytic decomposition of the perturbed system into two invariant 
blocks. By an argument similar to that for the geometrically regular case, we 
find that, at least for the easier nonglancing case, the associated block Gi^^ 
decomposes analytically as well. Thus, in the linearly splitting, nonglancing 
case, we may always arrange that the tangent system have unique (up to 
constants) symmetrizers. 

6.4 The case of double roots 

In this section, we assume that Ei is the direct sum of two dimensional 
subspaces Ei ^ associated to double linearly splitting characteristic roots: 

Assumption 6.15. Suppose that L is symmetric hyperbolic and nonchar- 
acteristic with respect to the boundary x = 0, and that its characteristic 
roots are either geometrically regular, totally nonglancing, or second order, 
nonglancing and linearly splitting transversally to a smooth manifold of codi- 
mension two. 

The following result is proved in Appendix D. 

Proposition 6.16. Suppose that Gi^k is a 2 x 2 block associated to a non- 
glancing linearly splitting eigenvalue, near (p, C) with 7 = 0. Then there 
exist smooth symmetrizers Si ^ near {p, (") for Gi k such that 



is a Friedrichs symmetrizer for L/^ j.. 

Moreover, the symmetrizer for L'^ ^ is unique, up to multiplication by a 
positive constant. 

In this case, the only block diagonal symmetrizers of L'^ are 



with positive constants Ofc. Moreover, given the a^, (|6.25|) defines a smooth 
symmetrizer of Gi such that 



Si,fc(p,C)A-i(p,C) 



(6.28) 





Therefore, Theorem 16. lUI takes the special form: 
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Theorem 6.17. Suppose that L satisfies Assumvtion \6. 1 ^ and M are bound- 
ary conditions which satisfy the uniform Lopatinski conditions. With no- 
tations as in Theorem \6.1(A there exists a smooth Kreiss symmetrizer for 
{G, M) if and only if for all {p, () with ( £ and j = 0, the reduced 
boundary conditions Mi {p, () are maximal dissipative for some Friedrichs 
symmetrizer o/Lp 

Moreover, we have by direct calculation the following useful fact, gener- 
alizing a corresponding observation of Majda-Osher in the case of the 2x2 
wave equation. For the proof, see Appendix D and Proposition 16.211 below. 

Proposition 6.18. If Gi itself is a 2 x 2 matrix, there is a smooth Kreiss 
symmetrizer for the reduced problem (Gi, Mi), if and only if it satisfies the 
uniform Lopatinski condition. 

This immediately implies the following 

Theorem 6.19 (Fourth Main Theorem). Suppose that Assumvtion \6.l3\ 

is satisfied and that for all C G S^ with 7 = 0, there is at most one real 
eigenvalue associated to a mode which is neither geometrically regular nor 
totally nonglancing. 

Then, for boundary conditions M there is a smooth Kreiss symmetrizer 
for (G, M) if and only if M satisfies the uniform Lopatinski condition. 

We conclude by briefly considering in the simple double-root case, under 
what conditions our structural hypotheses hold. 

Proposition 6.20. For a symmetrizable system in dimension d < 2, geo- 
metric regularity holds at any point 7^ 0. 

Proof. Restricting by homogeneity to the unit sphere, we obtain a symmetric 
matrix perturbation problem in a single parameter, for which eigenprojec- 
tions necessarily vary analytically |Katj . □ 

Proposition 6.21. In dimension d > 3, a variable multiplicity point of 
multiplicity two is generically nonglancing, nonregular, and linearly splitting 
of codimension two. Moreover, a linearly splitting point of multiplicity two is 
either of codimension two or else geometrically regular of codimension one. 

Proof. The reduced, symmetrizable system in dimension d — 1 = 2 gener- 
ically has noncharacteristic Ai, hence is nonglancing. Algebraic regularity 
can hold only if the reduced, frozen coefficient system has linearly varying 
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eigenvalues, in which case Ai and A2 commute, a degenerate case. More- 
over, if Ylj+i ^j^j equal eigenvalues for / (0, 0), then by symmetry it 
is a multiple al of the identity, which is three conditions (again taking into 
account symmetry to eliminate one entry) on the unknowns ^, a, generically 
determining a codimension two manifold M in 

The latter argument shows also that a linearly splitting point can be at 
most codimension two. Suppose that it is codimension one. Then, by an 
analytic change of coordinates taking the manifold M on which eigenvalues 
agree to the (Ci, • • • ,(,d-i) plane, and subtracting off the (analytic) average 
of the two eigenvalues, we obtain a matrix perturbation problem of form 
A(^) = ■^rf-D(^), where -0(0) is diagonalizable with distinct eigenvalues. The 
eigenvalues and associated eigenprojections of D are therefore analytic, as 
are the eigenvalues and eigenprojections for the original problem. □ 

Proposition 6.22. In any dimension, an algebraically regular variable mul- 
tiplicity point of multiplicity two is either geometrically regular or else has 
characteristics tangent to first order in ^, in which case nonglancing and 
totally nonglancing are equivalent. 

Proof. Denote the two characteristics as ti(^), T2(^). If 5(ri — T2)/d(,d) = 0, 
then one of the last two possibilities must occur. If 9(ri — T2)/dE,d) 7^ 0, 
on the other hand, then by the Implicit Function Theorem, there is a hy- 
persurface = ^Pi^') on which ti — T2 = 0. But symmetrizability im- 
plies semisimplicity, so that the reduced matrix A(^) obtained by projec- 
tion onto the associated two-dimensional total eigenspace of ri, T2 must be 
a multiple of the identity; likewise, dA/d^^d must be diagonalizable, with 
distinct eigenvalues. The matrix B{£^) := (^4 — ri/)(^) is thus of form 
{id - (t>{i'))D{^') + 0{\id - 'A(C')P)> analytic in ^, where D is diagonaliz- 
able with distinct eigenvalues, and it follows that the associated projectors 
are analytic. The final assertion follows by the same argument applied to 
each direction ^j. □ 

Collecting results, we have the following conclusion. 

Corollary 6.23. Assumption \6.15\ is satisfied always in dimension d = 
2. For multiplicity two crossings in dimension d = 3, assumption \6.1!^ is 
generically satisfied, and at algebraically regular points is always satisfied. 
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7 Application to MHD 



Finally, we discuss our main application, to the shock wave problem in MHD. 
As is well known, the Euler equations of MHD are symmetrizable hyperbolic 
for a thermodynamically stable equation of state, in particular for an ideal 
eras law; see ini iKawl iKShl iMZl IZ2] . However, they do not satisfy the 
condition of constant multiplicity, and so the standard Kreiss-Majda theory 
cannot be applied; see calculations. Appendix A. Moreover, the boundary 
condition that arises through linearization of the shock transmission problem 
is not dissipative, and so the older theory of dissipative boundary conditions 
cannot either be applied; see |BTlj . 

For this reason, and also the inherent complexity of the system, progress 
on shock stability for MHD has been so far somewhat limited. In particular, 
the only nonlinear stability result up to now has been obtained by Blokhin 
and Trakhinin |BTniBT4| using the "direct method" introduced by Blokhin 
|Blj for the study of gas dynamics, in the very special case that magnetic 
field H goes to zero: essentially, the fluid-dynamical limit. This consists 
of constructing a Lyapunov function in the form of certain "dissipative in- 
tegrals", from which linearized and nonlinear stability follow directly. The 
construction, which is rather intricate, depends crucially on symmetrizabil- 
ity of the system; as the authors describe it, it consists of augmenting the 
original system with various derivatives, chosen in such a way that the (ini- 
tially nondissipative) boundary condition becomes dissipative. At least as 
presented, it does not appear easily generalizable to other cases. 

On the other hand, with additional symmetry, either in the "transverse" 
case that the magnetic field be orthogonal to fluid velocity, or the "parallel" 
case that they are parallel, the Lopatinski condition has been explicitly 
evaluated even for large H, to yield physical stability conditions analo- 
gous to those of Majda |Maj| in the gas-dynamical case; see |GK1 IBTll 
IBT2| IBT3| IBT4) . More generally, the Lopatinski condition has been exten- 
sively studied numerically for arbitrary parameters, with interesting results 
|BTMH IBTM2j . For a detailed description of the current status of the 
theory, we refer to the excellent survey |BTlj . Up to now, however, the pre- 
cise relation between the Lopatinski condition and linearized and nonlinear 
stability was not known. 

Using the framework developed in this paper, we can now resolve this 
issue, at least for the simpler, isentropic case. We expect that similar con- 
siderations hold in the general case. 
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7.1 Equations 

The equations of isentropic magnetohydrodynamics (MHD) appear in basic 
form as 

{dtp + div(/9-u) = 
dt{pu) + div(pu*'u) +Vp + H X curlif = 
dtH + cui\{H X li) = 

(7.2) divF = 0, 

where p G M represents density, it G fluid velocity, p = p{p) £ R pressure, 
and H £ M!^ magnetic field. With H = 0, H7.1() reduces to the equations of 
isentropic fluid dynamics. 

Equations (|7.1|) may be put in conservative form using identity 

(7.3) H X curliJ = (l/2)div(|F|2/ - 2H^Hf'' + HdiwH 
together with constraint 1)7.2(1 to express the second equation as 

(7.4) dt{pu) + div(pu*'u) + Vp+ (l/2)div(|i7p/ - 2H*Hf' = 0. 

They may be put in symmetrizable (but no longer conservative) form by a 
further change, using identity 

(7.5) curl(F xu) = {divu)H + (n • V)H - {divH)u - {H ■ V)n 
together with constraint (|7.2() to express the third equation as 

(7.6) dtH + {di-vu)H + {u ■ V)H - {H • V)u = 0. 

7.2 Initial boundary value problem with constraint 

Constraint ()7.2() is preserved by the flow, hence may be regarded as a con- 
straint on the initial data. Thus, a sufficient condition for well-posedness 
under constraint (|7.2|) is that for some form of the equations, rewritten 
modulo (|7.2() . there hold well-posedness with respect to general initial data. 
Under the assumptions of the previous section, i.e. noncharacteristicity, 
symmetrizability, and the property that characteristics everywhere be ge- 
ometrically regular, totally nonglancing, or linearly splitting of order two, 
this is equivalent to satisfaction of the uniform Lopatinski condition for the 
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chosen representation of the equations. This is sufficient for our present 
needs, and seems to be the general approach followed in the literature. 

To obtain sharp (i.e., necessary) conditions, one should rather use the 
constraint together with the linearized equations to make a pseudodifferen- 
tial change of coordinates eliminating one variable, arriving at an uncon- 
strained system in one fewer variables. Equivalently, as described in |ChPj . 
one might by a pseudodifferential choice of augmented variables reduce the 
unconstrained second-order hyperbolic equation from which the system is 
derived to an unconstrained first-order system with modified boundary con- 
ditions. However, the resulting analysis would be special to the system 
considered, and not of the general form considered by Majda. 

In any case, our analysis here concerns sufficient conditions for stability, 
and so we shall follow the first, simpler approach and consider the equations: 



(7.7) 



' dtp + u ■ Vp + pdivu = 0, 

p{dtu + u ■ \7u) + Vp + H X cuilH = 0, 
^ dtH + uVH + {divu)H - H-Vu = 0. 



The Rankine-Hugoniot conditions are deduced form the the conservative 



form of the equations. Using notations (y, x) G 
as above, consider a shock front 



for spatial variables 



(7.8) 

with space time normal 

(7.9) (-cr,n), a = dt^p, n= {-dy^ip,...,-dy^_^ip,l). 

The jump conditions read, with Un ■= n ■ u and Hn = n ■ H, 

[p{un - a)] = 0, 



(7.10) 



[pu{un - a)]+n 



2' ' 



[HnH] = 0, 



[{Un - a)H] - [HnU] = 0, 
[Hn] = 0. 



The last jump condition comes from the constraint equation (|7.2|) . Appar- 
ently this system of 8 scalar equations is too large. However, projecting the 
third equation in the normal direction yields cr[i?„] = which is implied by 
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the last equation. This shows that H7.1U() is made of 7 independent equa- 
tions, as expected. Denoting by Utg and Htg the tangential part of u and 
H, that is their orthogonal projection on n-*-, ()7.1U|) is equivalent to 



(7.11) 



' [piUn ■ 
pu{Ur, 



a 



o-)l +n 



[HnH] = 0, 



[ [{Un - (T)Htg] - [HnUtg] = 0, [i/„] = 0. 

The constraint (|7.2() is preserved by the equations in the following sense. 

Proposition 7.1. Suppose that U = {p,u,H) is a smooth solution of (|7.7() 
for t £ [0, T] on both side of the front x = Lp{t, y) and satisfy the jump 
conditions 1)7. 11(1 . Assume in addition that the traces ofun — cr on both side 
never vanish on the front and at least one among the traces u~ — a and 
a — is positive. Then, if the initial value of H satisfies div//|(=o; then 
divH = for all t G [0,r]. 

Proof. The third equation in ()7.7|) implies that on both side of the front: 

dtH + curl(i7 X n) + {dwH)u = 0. 

Given a smooth function v on both side of the front, denote by v the distri- 
bution equal to v on both side of the front. The jump condition [i?„] = 
implies that divH = divH. Moreover, the jump conditions imply that 

dtH + cuvl{H X -u) = {dtH + curl(i? x u)y= -(dh^)M = -(div^)M. 
Therefore w := divH is a bounded and piecewise continuous solution of 



dtw + div{u'w) = 0, 



\t=o 



0. 



In particular, the jumps of u; = divH satisfy [(u„ — a)w] = 0, which could 
have been derived directly from the equations and the jumps conditions. 
Note that there are no general uniqueness theorem for equations dtw + 
div(au;) = with a G L°°. But here a is piecewise smooth ant the front 
is never characteristic for this equation since we assumed that u„ — a / 0. 
Moreover, the characteristics are impinging the front at least on one side. 
Therefore, integrating along characteristics, this implies that w = and 
thus that divH = 0. □ 
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7.3 The shock wave problem for MHD 

Thanks to Proposition 17.11 we consider the equations (|7.7jl on both side of 
a front ()7.8|) , plus the transmission conditions H7.11() . 

Fohowing the standard approach for hnearized stabihty analysis of fluid 
interfaces, (see |Maj| ) we carry the front location as an additional dependent 
variable, and make the change of (independent) variables 

X := X — ip{x' , t), 

to fix the shock location in the new variables at z = 0. In the new coordi- 
nates, the problem becomes 

2 

(7.12) dtU + J2^j(U)9yjU + MU,dip)di:U = 0, on ± x > 0, 
with 

2 

A3iU,dip) = A3{U)-J2i9yM^jiU) " (9*</^)Id. 
i=i 

The explicit form of the matrices Aj is given in Appendix A, where it is also 
verified that the system is hyperbolic symmetric in the sense of Friedrichs 
as soon as = dp /dp > 0. The unknowns are U = {p, u, H). The equations 
are supplemented with the transmission conditions (|7.11l) which read 

2 

(7.13) dtv[FQ{U)] +Y,dyMHU)\ = [F^m on x = 0. 
The stability analysis yields to consider the linearized equations 

2 

(7.14) dtU + Y,MU)dyfi + MU,d^)d3^U = f, on ± 5 > 0, 

2 

(7.15) dt^W^m+Y^dyA^jiU)] = [F!,{U,d^)U]+g on x = 0, 
with 

2 
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The Majda-Lopatinski determinant is computed as follows. One con- 
sider the linearized equations H7.14() ()7.15|) for U~ and constant and dip 
constant. This yields constant coefficients system, depending on parameters 
p = {U~ , , dip) to which we apply the analysis developed in the previous 
sections. Equations (|7.14|) (|7.15|) comprise a slightly nonstandard, constant- 
coefficient initial-boundary-value problem on the half-open space i > 0, in 
the "doubled" variables 

U+{z):=U+{z), U-{z):=U-{-z), 

the new feature being the appearance of front variable (p. As pointed out by 
Majda |Maj| , this may be handled within the standard framework, thanks 
to the cascaded form of the equations. 

For, assuming that the shock front is noncharacteristic, which follows 
from Lax' shock conditions, taking the Fourier transform in the transverse 
coordinates y and the Laplace transform in t, we obtain a system of ODE 
in variables U~^, U~ , 

(7.16) diU^ + G^{p,C)U^ = f^, on x^O, 
with as in 1)3.2(1 . and boundary condition , 

(7.17) 0X{p,C) = [FiU]+9, on x = 0, 
with 

2 

X{p,0 = {ir-j)[Fo{U)] + Y,vAFjm. 

i=i 

Noting that ip does not appear in the first equation, we may convert it to a 
separate, standard boundary- value problem by introducing a new, decoupled 
boundary condition of rank A^ — 1 = 6, obtained by taking the inner product 
of (|7.17|) with a basis of unit vectors in X{p, C)^. This reduces the boundary 
conditions to a system of the form 

(7.18) M{p,OU = g, ip = '^{p,C)U + gi, 

indicating how the boundary value problem analysis developed in the pre- 
vious sections applies to the present framework. 

More precisely, one introduces for 7 > 0, the spaces E±(p,C) of initial 
conditions corresponding to decaying solutions of (|7.16|) at ±00. They are 
generated by the generalized eigenspaces of G^{p, () located in {itlm /i < 0}. 
Lax' condition implies that 

(7.19) dim E~ (p, C) + dim E+ (p, C) = iV - 1 = 6. 
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The stabihty condition for the boundary condition 1)7. 17|) states that the ho- 
mogeneous equation with ^ = has no solution (ip, , U^) in C x x 
In accordance with the general Definition 14. 11 we obtain Majda's Lopatinski 
determinant for the full shock problem, 

(7.20) D{p,C) =det(^F^-]E-, Ff^+E+ , Cx) , 
or 

(7.21) D(p,C) =det(rf,...,r-_i, rp+i,...,rf, X) , 
where rj and are bases for E~ and E"*". 

7.4 Applications of Theorem 15.61 

Following the analysis presented in the previous subsection, we consider the 
linearized equations 1)7.14(1 (|7.15() around a planar shock. Without loss of 
generality, we can assume that the spatial normal to the front is n = (0, 0, 1) 
so that the normal matrix is = — aid. The seven eigenvalues of 
Ml], = m^i + 112^2 + are 

Ao = r]iui + r]2U2 + (us - o-)C 

f7 22) A±i = Ao±c.(H2+a^/^ 

I A±2 = Ao ± imHi + mH2 + ^H^)/^ 

^A±3 = Ao±c/(|77|' + e')'/' 

with 

cj := i (^c2 + h^) + y/(c2 - /i2)2 + 452^2^ 
c2 := i (^c^ + h^) - V(c2 - /i2)2 + 452^2^ 

with c2 = p'{p) the sound speed, h'^ = \H\'^/p, b'^ = |(r/,C) x H\'^/p and 
(^7)0 = (^) 0/1(^5 01 (see Appendix A for the computations). The first 
eigenvalue corresponds to the transport of the constraint. It can be decou- 
pled from the system : there is a smooth one dimensional subspace. Eg such 
that A{ri,^) = Xq on this space and Eq- is stable for A{r],^). The other 
eigenvalues are in general simple. 

The boundary {x = 0} is noncharacteristic if and only if 

(7.23) us - cT ^ {0, ±Hs/^, ±csin), ±c/(n)} 

where Cs{n) and Cf{s) are the slow and fast speed computed in the normal 
direction n = (0,0, 1). The next lemma is proved in Appendix A. 
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Lemma 7.2. Assume that < \H\'^ ^ pc? . Consider {r],^,) in the unit 
sphere 5^. 

i) When (jy,^) • / and (r/, ^) x H ^ 0, the eigenvalues are simple. 

ii) On the manifold {r],^^) ■ H = the eigenvalues X±s are simple and 
the multiple eigenvalue Xq = \±i = \±2 is geometrically regular. 

Hi) On the manifold (r/,^) x H = 0, Xq is simple. When \H\'^ < pc^ 
[resp. \H\'^ > pc^ ], [resp. \±i are simple, the other eigenvalues A+2 
and A_2 o^'e double, algebraically regular, and totally nonglancing provided 
that — a ^ zb^Ts/y^. 

This shows that the assumptions of Theorem 15.61 are satisfied as soon as 
the hyperbohcity condition (? := dp /dp > holds, and H ^ with length 
\H\ ^ cy^ and condition (|7.23|) is satisfied. 

Corollary 7.3. Under the generically satisfied conditions < 7^ 
/9^(c^)^, the uniform Lopatinski condition implies linearized and nonlinear 
stability of noncharacteristic Lax-type MHD shocks. 

Proof. Lemma 17.21 and Theorem 15.61 imply the existence of smooth Kreiss 
symmetrizers. The stability follows through the standard shock stability 
framework of Maj} IMe2j (see the useful Remark 15. 9|) . □ 



7.5 The H ^0 limit 

Using Corollary 17.31 we easily recover the result of Blokhin and Trakhinin 
in the isentropic case, by a simple perturbation argument. When H = 0, 
the system (|7.7|) reduces to isentropic Euler's equations and ()7.11|) to the 
corresponding Rankine Hugoniot condition. Denote by Do{po, the Majda- 
Lopatinksi determinant of this problem, where the parameters are now pQ = 
{Uq , Uq , dip) with Uo = {p, u). Note that because we have everywhere either 
geometric regularity or total nonglancing, the stable subspace E_ , and thus 
the Lopatinski determinants have continuous extensions to C G (see 
Theorem 15. 6() . 

When H = 0, the eigenvalues are 

, . Ao = A±i = A±2 = Viui + + ^(^3 - cr), 

A±3 = Ao±c|(r?,0|. 
The boundary is noncharacteristic when 
(7.25) U3-a {-c,0,+c}. 
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Lemma 7.4. Suppose that U_ = (p, n, 0) and a satisfy ()A.26|) holds. Then 
the multiple eigenvalue Aq is totally nonglancing at (U_, rj, ^) for all (r/, ^ 
and the eigenvalues A±3 are simple. 

The proof is given in the Appendix A. By Theorem 15. 61 this imphes that 
the negative spaces depend continuously of H for H close to 0. Therefore, 
the Lopatinski determinant D{p,Q is a continuous function of G S_^_ and 
p = {U-,U+,df), with = {p^,u^,H^), when C/± satisfy ((7^^ and 
are small. Therefore: 

Corollary 7.5. Suppose that po is a Lax shock for the insentropic Euler's 
equation. Then, as {H^,H^) 0, there holds D{p,Q Dq{pq^C,), uni- 
formly for £ S^, with 7 > 0. 

Corollary 7.6. In the H ^ limit, Lax-type MHD shocks approaching 
a noncharacteristic limiting fluid- dynamical shock satisfy A ssumption \6.L% 
so that the Lopatinski condition reduces to 1)4. lOj) .' moreover, they satisfy 
the Lopatinski condition if and only if it is satisfied by the limiting fluid- 
dynamical shock (i.e., satisfies the physical stability conditions of Erpenbeck- 
Majda \Eri \Majli , in which case the MHD shocks are linearly and nonlinearly 
stable. In particular, for an ideal gas equation of state, they are always stable 
in the H ^ limit. 

Likewise, using Corollary [721 we may immediately convert the results of 
[(irKj . |BTlj -jBT4 . etc. on satisfaction of the Lopatinski condition to full 
results of linearized and nonlinear stability. Adore generally, Corollary 17.31 
gives a useful framework for the systematic study of AIHD shock stability 
and other variable-multiplicity PDE problems through the investigation of 
the linear-algebraic Lopatinski condition. 

Remark 7.7. Corollarv 17.61 guarantees stability for sufficiently small but 
nonzero values H^. But since the system is symmetric hyperbolic for all H, 
a direct application of Theorem 15.61 provides us with Kreiss symmetrizers 
which are smooth in p = {U~ , U~^,d(p) for \II^\ small. Therefore, we have 
maximal stability estimates for the linearized equations which are uniform 
in for small. This implies uniform nonlinear stability, since the time 
of existence depends only on the estimates (see Maj| , |A/[e2) ) . This is likely 



to imply the continuity of the shock solution as -ff — > 0. 

Remark 7.8. Calculations as in |Me2llZ2j show that Lax-type MHD shocks 
are stable in the small-amplitude limit [U] — > 0, again, for H± ^ 0. 
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A Appendix A. The symbolic structure of MHD 
A.l Multiple eigenvalues 

The first order term of the linearized equations of ()7.7() about {u, H) is 



(A.l) 



' Dtp + pdivu 
Dtii + p^^c^Vp + p^^H X curliir 



, DtH + (divn)F -H -Vii 

with Dt = dt + u ■ V and = dp /dp. This system is hyperbohc symmetric 
with symmetrizer S = block-diag(c^, /3ld, Id). The associated symbol is 

(A.2) 



' TP + p{i ■ ii) 

fii + p'^c^pi + p'^H X {^x H) 
[fH + {C- u)H - {H ■ e)n 



with f = T + u ■ We use here the notation ^ = (■^1,^2,^3) foi^ the spatial 
frequencies and 



U\\ 



C • u , U± = U — Ull^ = X X u) 



We write HA.2|) in the general form rid + A{U, ^) with parameters U 
{p,u,H). The eigenvalue equation A{U,S,)U = XU reads 



c^p + H^-H^, 



(A.3) 



' Xp = pii 
pXii 

pXu± = —H\\H±, 
XH± = liiii^j, — H\\u^, 
XH\\ = 0, 



with A = A — (n^). The last condition decouples. On the space 
(A.4) Eo(e) = {p = 0, n = 0, i7x = 0}, 

A is equal to Xq := u ■ ^. From now on we work on Ej^ ~ {^|| =0} which is 
invariant by A{p,^). 
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Consider v = H/^/p, v = H / ^ and a = p/ p. The characteristic system 
reads: 



(A.5) 



' A(T = till 
\u\\ = (?a + v±- v±_ 

, \v± = U\\V± — V\\U± 

Take a basis of such that v±_ = (6, 0) and let a = v\y In such a basis, the 
matrix of the system reads 



(A.6) 



A 


-1 
















A 








-h 











A 





a 














A 





a 





-h 


a 





A 














a 





A 



The characteristic roots satisfy 

(P - a2)((A2 - a2)(A2 - - AV) = . 



(A.7) 




- a 


Thus, either 






(A.8) 


A2 = 




(A.9) 


A2 = 




(A.IO) 


A2 = 





C~ + h^) - V(c2 - /l2)2 + 462c2 



2 
2 



with h? = a'^ + 11^ = \H\'^/p. 

With P(X) = (X-a2)(X-c2)-62x, {P < 0} = [cl,c)] and P(X) < 
for X G [min(a^, c^), max(a^, c^)]. Thus, 



> max(a^,c^) > a 



(A.ll) 

(A. 12) cl<mm{a^,c^) <a^ 

1. The case t'^ 7^ i.e. w = ^ x ?; 7^ 0. Thus, the basis such that HA.6|) 
holds is smooth in In this basis, w = {d,h) , h = \v ±\ > 0. 
1.1 The spaces 

E±(|) = {(^ = 0, n|| =0, {;_L G C(| x v), n_L = ^F-^^} 

are invariant for A and 

(A. 13) A = ±a on E± . 



60 



1.2 In (E+ © E-)-*-, which is invariant, the matrix of A is 
(A.14) Ao := 



10 

-6 

-a 

-6 a 



Since P{(?) = -b'^c^ < 0, there holds < < cj. 

1.2.1 ) Suppose that a 7^ 0. Then, ^(a^) = -a^c^ < and 
c1 < < cj. Thus, all the eigenvalues are simple. Moreover, c^cj = a^c^ 



and > 0. The space 



Attii, ttii 



A2 



A 



is an eigenspace associated to the eigenvalue A when A = itc/ and A = ztc^. 
Here iii and vi denote the first component of u± and v± respectively in the 
basis {v±,w). 

1.2.1 ) Suppose that a is close to 0. Since cj > > 0, the spaces 
¥±cf are still eigenspaces associated to the eigenvalues A = itc/. 
By direct computations: 



;2 + /i2 



Therefore, 



> as a — s- ., . 



Therefore, consider Cg = j^Cs is an analytic function of a (and b ^ 0) near 
a = and 

F±,,(a,6) =F±e, 

are analytic determinations of Eigenspaces, associated to the eigenvalues 
zbcg. Moreover, the values at a = are 



F±,.(0,6) = {(7 = 0, ny =0, u± = 



and F+,s n F 



{0}, thus we still have an analytic diagonalization of Ao- 



61 



2. Suppose now that b is close to zero. At 6 = 0, the eigenvalues of A 
are ±c (simple) and ±h (double). Assume that ^ h^. Note that when 
6 = 0, then \a\ = h and 

when cP > h'^ : Cf = c, Cg = h, 
when < : Cf = h, Cg = c. 

2.1 The eigenvalues close to ±c remain simple. 

2.2 We look for the eigenvalues close to h. The characteristic equation 
implies that 



(A.15) 



{c^a = A^ii — v± ■ v± 
(A^ — c^)u\\ = Xv± ■ v± 

Eliminating u\\, we are left with the 4x4 system in ^-^ x 

Anj_ = —av± 



(A.16) 
Thus, 



-au± + - 



A 



(P - a^)v^ 



A2-c' 
A2 



:{v± (g) v±)v±. 



A2 



{v± (g) v±)v±. 



Recall that \v±\ = 6 is small. We recover 4 smooth eigenvalues 

(A.17) ±a, ± + 0{P) = ±(a + 0(6^)) . 

(remember that a = ±/i + 0(6^). However, the eigenspaces are not smooth 
in V, since they are Rv± and x v± and have no limit aX v± ^ 0. 

Summing up, we have proved the following. 
Lemma A.l. Assume that = dp/ dp < 0. The eigenvalues of A(U,^) are 



(A.18) 



Ao = ^ • n 

A±i = Ao±c,(C)|C| 
A±2 = Ao ± (^ H)/^ 
U±3 = Ao±c/(|)|^| 
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with ^ = C/ICI o.'^'d 



(A.19) 




(A.20) 



w/iere /i^ = = |l x 

Lemma A. 2. Assume that < 7^ pc^ where = dp/dp > 0. 

When ■ H ^ and x H ^ 0, the eigenvalues of A{U,^) are 

simple. 

a) On the manifold ■ H = 0, ^ 7^ 0, the eigenvalues A±3 are simple 
and the multiple eigenvalue Aq = A±i = A±2 is geometrically regular. 

Hi) On the manifold ^ x H = 0, ^ ^ 0, Xq is simple. When \H\'^ < pc^ 
[resp. \H\'^ > pc^ ], A±3 [resp. \±i] are simple; A4-2 7^ A_2 a'^e double, equal 
to \±i [resp. A±3 ] depending on the sign of^-H , algebraically regular, but 
not geometrically regular. 

A. 2 Glancing conditions for MHD 

Write the space variables (y, x), y G M^, x G M and consider the planar front 
{x = at}. As indicated in Section 7, we perform the change of variables 
X = X — at, obtaining the linearized system H7.14() with reads 




with U = {p,u, H), u = {ui,U2,U3 — a). This shows that the eigenvalues of 
the symbol are 



where {r],^^) denote the tangential Fourier frequencies dual to {y,tildex). 
Therefore, Lemma lA. II implies that the eigenvalues are given by H7.22() . 
In particular, the eigenvalues of the boundary matrix are 

(A. 23) ns-fj, U3 - a ± Hs/^, U3-a±Cf{n), U3-a±Cs{n). 

with n = (0,0,1). The boundary is noncharacteristic if they are different 
from zero, that we assume from now on. The next result finishes the proof 
of Lemma l7. 21 

Lemma A. 3. Suppose that < \H\'^ ^ p(? and the front in noncharacter- 
istic. On the manifold (rj, x H = 0, (rj, ^) 7^ 0, the double eigenvalues A+2 
and A_2 o,i"s totally nonglancing. 




(A.21) 



dt - adi + A{U, dy, 9,) =dt + A{U, dy, d,) 



(A.22) 
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Proof. Consider 



(A. 24) A2 = mifJ-i +vi) + 772(^2 + V2) + C{u3 + V3-a), 

with V = H/^/p. We have shown at (|A.17|) that the eigenvalue is alge- 
braically regular, and the nearby eigenvalue is 

A2 + 0(|C X 

It is tangent to A2 to second order. Therefore, the nonglancing condition 
reduces to 

— - = u3+V3-a j^O, 
that is automatically satisfied when the front is noncharacteristic. □ 
We now pass to the proof of Lemma 17.41 When H = 0, the eigenvalues 

are 

. Ao = A±i = A±2 = Viui + mu2 + ^(^3 - cr), 

A±3 = Ao±c|(r/,0|. 

The boundary is noncharacteristic when 
(A.26) U3-a ^ {-c,0,+c}. 

Lemma A. 4. Suppose thatU_= (p, u, 0) and a satisfy ()A.26|) holds. Then 
the multiple eigenvalue Aq is totally nonglancing at (U_, r], ^) for all (r], (,) ^ 0. 

Proof. Note that Aq has constant multiplicity 5 in (77, ^) for H = 0, and 
is linear in Therefore, the tangent system at (f/, r/, ^) is AqH in 

dimension 5. In particular, with notations as in Section 3, the boundary 
matrix of this system is -A3 = {u3 — cr)Id, implying that the eigenvalue is 
totally nonglancing. □ 

B Appendix B. Maxwell's equations in a bi-axial 
crystal 

We give here an example where geometric regularity fails. The consequence 
is that the tangent operator is no longer a transport operator along a single 
direction, but rather a system of hyperbolic equations of wave equation 
type. This example permits the description of the phenomenon of conical 
refraction (cf [HlTil F^iyi IRW] 'l 
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Maxwell's equations in the absence of exterior charge may be written as 



(B.l) 



dtB + curlE = , divB = , 
dtD - cmlH = , divD = , 



with 

(B.2) D = SE, B = nH 

where £ and /x are 3x3 positive definite matrices. In the case of a bi-axial 
crystal, fi is scalar and £ has three distinct eigenvalues. We change variables 
so that /X = 1 and choose coordinate axes so that £ is diagonal: 

/ a, 

(B.3) £-^ = I ^2 

V 

with ai > 02 > 03. Ignoring the divergence conditions, the characteristic 
equation and the polarization conditions are obtained as solutions of system 



E J ' \tE - £-^{^ X B) 

For 7/ 7^ 0, r = is a double eigenvalue, with eigenspace generated by (77, 0) 
and (0, ry). These modes are incompatible with the divergence conditions. 
The nonzero eigenvalues are given as solutions of 

E = £-\^xB), {t^ + n{0£-^n{^))B = 

where is the symbol of the operator curl. We have 

(-0:2^3 "366 "266 
"266 -aiCl - «2^? 

det(r2 + A{0) = t^{t' - ^{iV + \e m) 

with 

f ^{i) = (ai + a2)e| + (a2 + ag)?! + (a3 + "i)^! 
1 ^(0 = "102^3 + "203^ + OL^aiil ■ 
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The equation det(T^ + A{^) = yields r = and an equation of second 
order in , of which the discriminant is 

with 

P = (ai- 02)^1 + ("3 - "2)^? + (03 - "1)^2 
Q = 4(ai-a2)(ai-a3)e|ei>0. 
The root is double if and only if 

(B.5) 6 = 0, aiCi + «3^? = MCl + fs) = ■ 

Let (r, ^) be a solution of ()B.5|1 . For 6 such that 03 cos^ 9+ai sin^ 9 = a2, 
we have 

(B.6) i = Xbi, T = ±X^. 

where we have introduced the basis of 



hi 



In HB.6|) . suppose now that r = Xy/012 , the other case being similar. The 

kernel of L(r, ^) is of dimension two, generated by 

(B.7) 

and 5j := ajj^/ai- As L{t,S,) is self-adjoint, (^2,^3) also form a basis of 
ker L*(t,^). The tangent system L/ at (t, ^) is obtained by multiplying on 
the left by a basis for the left kernel and on the right by a basis for the right 
kernel of L(r, 1^), see Remark 12.71 We obtain a 2 x 2 system with symbol 

/g / T- (,fi(53Cos6' + ^3(5isin6') 6 (5i - (^3) cos 6* sin 6* 

V .^2(<^i — (^3) cos6'sin6' r — (,^i(52 cos + .^3(52 sin ( 

Compare to equations 1)2. 25() and (|2.26|) obtained in the semisimple and 
geometrically regular cases. By a change of variables, we arrive at the form 

(B.9) (^~^^' 
This shows that the eigenvalue is linearly splitting. 





( -5isin0 






with 62 = 




1 — 1 






I 82, cos 9 ) 
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C Appendix C. The block structure condition 

We consider the system 

(C.l) d, + iG{p,0 

With ^ = (r — ^7, 77). Assume that 

Assumption C.l. The characteristic polynomial in ^, A{p, (, ^) = det (^Id+ 
G{p, Q), has real coefficients when 7 = and has no real roots when 7 > 0. 

Definition C.2. G has the block structure property near {p, () if there exists 
a smooth matrix V on a neighborhood of that point such that V^^GV = 
diag(Gfc) is block diagonal, with blocks Gk of size X Vk, holomorphic in 
T — ij, having one the following properties: 

i) the spectrum ofG^ipX) contained in {Im/x / 0}. 
a) Uk = 1; Gk{p, C) is real when when 7 = 0, and d-yGk{p, C) / 0, 
iv) > 1; Gk{p,C) has real coefficients when 7 = 0, there is fi^^ & M. 
such that 



(C.2) Gfe(p,0=/xjd + 



and the lower left hand comer of d^Qkip,C) does not vanish. 

Theorem C.3. The block structure condition is always satisfied when 7 > 0. 
When 7 = 0, it holds if and only if for all real root of A(p, ^) = 0, 

i) there are smooth functions Xj{p,r],$,), analytic in real when ^ is 
real, near (j>,r],$,) such that 

i 

(C.3) A{p, C, = eip, C, 0ll{r-i7 + Xjip, rj, 0) 

j=i 

with e{p, C, 7^ o'^^ ^j{Pj + r = 0, 

a) there are smooth vectors ej{p, ry, ^) near {p, rj, Q, analytic in ^, linearly 
independent, such that 

(C.4) (CId + G(^)(p,r?,e))e,(p,r/,e) = 

where G^^^ is the matrix G evaluated at t — ij = — Aj(p, r/, 



1 

■•• 

■•. ■•• 1 

... 
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Proof. The first statement is clear from Assumption IC.ll Moreover, one 
can always make a first block diagonal reduction with blocks associated 
to distinct eigenvalues of G{p,(^). Non real roots always satisfy condition 
i) of Definition IC.2I So is is sufficient to consider separately each block, 
associated to one real eigenvalue. Thus, it is sufficient to prove the theorem 
when 

(C.5) eid + G(p,C) = 0, eeM. 

a) Suppose that the block structure condition holds. Then 

(C.6) A{p,C,C) = llAk{p,C,0, Afc(p,C,0 :=det(eid + Gfc(p,C)). 

The form ()C.2|) of Gk, implies that 

det (e + Gp,r-i7,r/)) = {-l^'^-^c^ + 0{j^) 

where c is the lower left hand corner of djQk{p,C)- Since c 7^ and Qk is 
holomorphic in t — i'j, this implies that 

drAk{p,C,i)y^O- 

Since is real when and t — i-f are real, the implicit function theorem 
implies that there is a smooth function Afc(p, r/, ^), analytic in ^ such that 

(C.7) Ak{p,C,0=e{pX,0{r-il + Xk{p,r],0 

on a neighborhood of (p, (, ^), with e{p, CO / 0- 

Decompose the matrix V into (^i, . . . , 14) where T4 is of size N xu^ and 
denote by V^''\p,r],^) the matrix Vk evaluated at r — 27 = —Xk{p,r],^). At 
the base point 

(C.8) vPip,ri,i) = Vkip,C)- 

If i^k = 1, then the coefficient £,+QkiPj C) vanishes at when t — i"f = —\k 
and Tfc = vj^^'^ is a smooth element in the kernel of + G^'^^p, rj, 

If z^fc > I7 consider next the first z^fc — 1 rows of the eigenvector equation 
(^Id+Gjt)n = 0. Denoting by vi the first component and by v' the remaining 
ffc — 1, these equations read 

(C.9) Q{p,C.€)v' = viKpX.O 
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with C)0 = Id s-iid C)0 = 0- Therefore, one can solve v' = 
viQ^^h. When t — = — Afc(p, rj, the determinant of ^Id + G^^\p, rj, 
vanishes. This shows that the kernel of ^Id + G^^^ {p, rj, ^) has dimension one 
and is spanned by the smooth vector = '^{l,Q~^h) where the functions 
are evaluated at r — 17 = — Afc(p, 77,^). This implies that := V^^^ek is 
a smooth vector in the kernel of ^Id + G^^^. By ()C.8|) . the vectors are 
linearly independent at {p, r], and thus remain independent in the vicinity 
of that point. 

b) Conversely, suppose that (|(7.4|) and ()('.5|) are satisfied. By ()('.5|) . 
there is only one real eigenvalue — ^ at (p, C)- For all j, there is an integer 
Uj > 1 such that 



(C.IO) d^Xj = ... = dl' ^Xj = 0, dl'Xj / at (p,r/,0. 



Repeating the analysis in |Me3j . this implies that 

(C.ll) r-i7 + Aj(p,r/,e) = ej(p,C,e)Ai(p,C,e) 

with ej{p,C,0 / and Aj{pX,S,) ^ monic polynomial in ^, of degree vj, 
with real coefficients when 7 = 0, such that 

(C.12) A,(p,c,e) = (c-ir. 

Suppose that Uj = 1. Then, Aj = ^ — fJ,j{p,C) with fij smooth and 
^jipX) = ^jiViVi ^''j) is a smooth eigenvector of G{pX) with eigenvalue ^j. 
Comparing with ()C.11|) . implies that d^Hj{p,Q 7^ 0. Thus, the one dimen- 
sional space Cej provides us with a one dimensional block which satisfies 
property ii) of Definition IC. 21 

Suppose next that Vj > 1. Differentiate the equation HC.4I) with respect 
to ^, at Since G'^^\p,ri,^) = G{p,C) + 0((x - CY^), this implies 

that 

(C.13) e,., = i^(a|e,)(p,r/,0 

satisfy 

(C.14) (e + G(p,C))e,,o = 0, 

{§_ + G(p, 0)%, = fc-i) fc = 1, . . . , i/j - 1. 

The definition of the ej^k is extended to {p, C) in a neighborhood of (p, (), as 
in [Me3] : 
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where r > is such that Aj(p, ^) has no root on the circle {|^— ^| = r}. By 
HC.12() . the two definitions agree at {p, (), so that ej^kip, C) = §.j,k- Repeating 
the analysis of |Me3j , one shows that the space Kj (p, C) spanned by the 
Bj^kip, C) is invariant by G{p, (). It is smooth and of dimension Uj. 

Moreover, using ()('.14|1 and the independence of the ej{p,C), one shows 
that Cj^fc, for A: G {0, ... , Uj — 1} and j G {1, . . . , j} are linearly independent. 
Therefore, this provides a smooth decomposition of into invariant sub- 
spaces of G, thus a block decomposition of G. Moreover the characteristic 
polynomial of Qj is Aj. The property ()(L2|) is stated in (|(].14|) . 

The proof goes on as in |Me3j . For each block, one can perform an 
additional change of bases such that 



Qk = Qk{p,0 + 



qi{z) 
qv{z) 








(see (Rail). Thus, 

V 



v-l 



1=1 



Since Aj has real coefficients when 7 = 0, this implies that the qi{z) and 
therefore Qkip, C) is real when 7 = 0. In addition, by (jC.lip . 

Therefore, Q{z) satisfies the property iv) of Definition IC . 21 and the proof of 
Theorem IC.3l is complete. □ 



D Appendix D. 2 x 2 sytems 

Consider the 2x2 complex system in space dimension 2: 
(D.l) rid + ^A{p) + r]B{p) 

Assumption D.l. The system is strictly hyperbolic in the direction dt and 
noncharacteristic and nonhyperbolic in the direction dx. 



70 



D.l Reductions 

In particular, the eigenvalues of A are real and distinct. Thus, there is 
a smooth change of basis such that A is diagonal, with nonvanishing real 
diagonal coefficients ai ^ a2- One can perform changes of variables which 
preserve the directions dt and dx: 

(D.2) t' = T + ari, ^' = ^ + /?r/, t]' = rj. 

This transform HD.lj) to 

(D.3) T'ld + i'A{p)+ri'B'{p). 

With appropriate choices of a{p) and one can cancel the the diagonal 
terms in B': 



(D.4) B' 



b 
c 



Strict hyperbolicity implies that be is real and positive. By a change of basis, 
one can assume that 6 = c is real, so that the system reads 

(D.5) r'U + d""' ' 



' 0, 

The assumption implies that aia2 < 0. With a = (—01/02)^^^ > and 

(D.6) r" = r', e" = -e', v" = b^, 

ai 

we obtain the reduced form 

It remains only one parameter a = a{p) > 0. 
D.2 Negative spaces 

Dropping the " we consider ()D.7() . and more precisely 



(D.8) eid + G(p,r,ry) :=ad + 

The eigenvalue equation for G is 



r/a r]/a 
—ar] —ar 



(D.9) + aTn--Tfi = T'^ -r]'^. 



a 
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If ^ + ar / 0, the corresponding eigenspace is 



(D.IO) E = Cr, r = f ,/ 

For Imr < 0, there is one eigenvalue /i in Im/i < 0. We denote by E_(r, ry) 
and r(r, ry) the associated eigenspaces and eigenvectors. 

Let z = r]/{fj. + ar). Then, for z ^ 0, (|D.9|) is equivalent to 

(D.ll) (a + -)T = ri(z + -). 

a' ^ z' 

In particular, 

(D.12) (a + i)ImT = T/Imzfl - -i^ 



a' ^ ""i^ 



When Imr < and 77/0, Im(/i + ar) < and Imryz > 0, and (|D.12|) 
implies that \z\ < 1. 

Conversely, for \z\ < 1, Imz 7^ 0, we can choose = ±1 such that 
rjlmz > 0, implying that r defined by HD.11|) satisfies Imr < 0. Then 
fj, = —ar + rj/z is a root of (|D.9() and 



r/Imz T/Imzl + a^lzP 
Im/i = -almr — -r- = — ; — < 0. 

Lz r z r 1 + 



This shows that the mapping (r, r/) 1-^ z sends {Imr < 0;r/ 7^ 0} onto 
{\z\ < l;Imz 7^ 0}. By continuity argument, or by direct computation, this 
implies the following. 

Lemma D.2. The union of the complex lines E_(r, r/) for Imr < and 

{T,r]) ^ is the cone T C of vectors u = *(Mi,n2) such that a\u2\ < 

Corollary D.3. The boundary condition Mu = ui — cu2 = satisfies the 
Lopatinski condition, if and only if a\c\ < 1. 



D.3 Symmetrizers 

In the form HD.7|) . or HD.5|) there is a unique symmetrizer (up to a scalar 
factor) which is the identity matrix. Tracing back to the original system 
(|D.1|) . yields a unique (up to scalars) symmetrizer S{p) = V*{p)V{p)S_ where 
V is the change of basis. The symmetrizer for G is S = SA. 

Proposition D.4. The boundary condition Mu = satisfies the Lopatinski 
condition if and only if it is maximal dissipative. 
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Proof. This is invariant by the change of basis. For systems HD.7|) . maximal 
dissipative and uniform Lopatinski conditions are such that kerM = {ui = 
CU2}- Since 

= a|ni|^ W2\'^, 

a 

the boundary condition is maximal dissipative if and only if this form is 
definite negative on kerM, that is if and only if a^|cp < 1. □ 

D.4 2x2 linearly splitting blocks 

We are given a 2 x 2 matrix G{q, r, t/, 7), = (r, rj, 7) in a neighborhood of 
in M^, and q a set of parameters in a neighborhood of q. with the property 
that 

(D.13) G{q,C)=(p{q)ld + G{q,C), G{q,0)=0. 

Assumption D.5. i) The characteristic polynomial in ^, A{q, ^) = det ('^Id+ 
G{q, C), has real coefficients when 7 = 0. 

a) The first order expansion in ( of G at q is 

(D.14) G'(C) = (f-i7)Go + ^Gi 

where detGg < and rid + Gj^^{C + ^G^) is strictly hyperbolic. 

Proposition D.6. There is a smooth self adjoint matrix C) on a neigh- 
borhood of {q, 0) such that 

i) — Im(SG) = ^yE, with E{q,0) positive definite, 

a) with S = 0), SG q and T,G i are self adjoint. 

Proof, a) We can assume that Gq is diagonal. Introduce the notations 



The assumption implies that 

«o ~ 9ra{q, 0) and dg = drb{q, 0) are real and agdg < 0, 
(D.16) ~ ~ 

bi = drib{q,0), £1 = dric{q,0), b^Ci is real and b^Ci < 0. 

By assumption, there holds 

(D.17) a + deR, ad-bceR when 7 = 0. 
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We look for a symmetrizer 

(D.18) ^"(^ 5)' a,5eR,l3eC. 

b) We show that one can choose such that 
(D.19) Im(SG)=0 when 7 = 0. 

In this part, we assume that 7 = 0. ()D.19|) is equivalent to the conditions 

(D.20) lm{aa + pc) = 0, 

(D.21) lm(pb + 6d) = 0, 

(D.22) {a-d)p = ab-6c. 

Because of ()D.17|) . a := a—d = Re (a—d) is real and by HD.16|) . drcr{q, 0) 7^ 0. 
Therefore, we can take {q, a, r/) as local coordinates near {q, 0, 0), and we first 
solve the equation 

(D.23) ab-Sc = 0, on cr = 7 = 0. 

By (ITTTbI) . there holds 

(D.24) V=o,7=o = Vh + 0(7?2), C|,=o,7=o = VQi + 0(r?^) 

Moreover, when a = 0, ad = |ap and by HD.17|) . be is real. Therefore, 
factoring out ry in (|D.2,3|) . we see that there is a smooth local solution a, 6, 
such that 

(D.25) a((?,C)GM, <5(g,C)eIK, a{q,0) = \ci\^ 6iq,0)=hc,. 

By ()D.23|) , one can factor out a in ab — 6c restricted to 7 = and there is a 
smooth function (3 such that HD.22|) holds when 7 = 0. In particular, since 
b and c 

(D.26) /3(g,0) = 0. 

For 7 = and o" 7^ 0, 

a ct 
Im (oa + /3c) = — Im (aa + be) = — Im (6c — ad) 
a a 

vanishes by HD.17|) . Thus, HD.20|) and similarly, (|D.21I) are satisfied on 7 = 0. 
Therefore, we have constructed a smooth solution of (|D.19|) . 
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Morever, differentiating this equation with respect to r and ry, imphes 
that 



(D.27) Im (SGo) = 0, Im fSGj ) = 0. 

The first equahty, ^PJK^ and (|Dl6|) . imply that P{q,0) = 0. Thus, with 
HD.25|) . we see that 

(D.28) - = (o 5) ' a5<0. 

c) Changing E to — E if necessary, we can assume that ociq > 0, implying 
that Mq > 0, hence that EGq is definite positive. 

The identity HD.19|) implies that Im (EG) = —'jE with E depending 
smoothly on {q, () in a neighborhood of {q, 0). Moreover, at (7 = g, r = = 0, 
there holds 

Im (EG) = Im (-^tEGq) + 0(7^) 
showing that E{q, 0) = EG n is definite positive. □ 

Theorem D.7. Suppose that Assumption \D. ,5l is satisfied and that M{q,C) 
is a 1 X 2 matrix such that the equation 

(D.29) d,u + iG{qX)u = f, M{qX)u\,=o = 9 

satisfies the uniform Lopatinski condition near {q,0). Then, there is a 
smooth symmetrizer E((/, C) near (q, 0) and there are constants c > and 
C , such that 

(D.30) -lm(EG) > C7ld, when 7 > 0, 

(D.31) -E + CM*M > cid. 

Proof. The symmetrizer E has been constructed in Proposition ID. 61 To 
prove that HD.31|) holds on a neighborhood of ((;, 0), it is sufficient to prove 
that it holds at {q, 0), thus it is sufficient to check that — E is definite positive 
(or E definite negative) on ker M, where M = M(g, 0). 

In addition, since EGg is self adjoint and definite positive, and EG ^ 
is self adjoint, E is a symmetrizer for G' = (t — i7)Gg + r]Gi. Thus, by 
Proposition ID. 41 it is sufficient to prove that the boundary condition M 
satisfies the uniform Lopatinski condition for the operator dx + iG'. 

This follows from the analysis in Section 4. We have shown that the 
negative space for G{q, C) with 7 > 0, is the negative space for a matrix 
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G{q,p,C) with p = Id and C = C/KI- Moreover, ^(^.O,^ = G'(C). The 
Lopatinski condition imphes that 

(D.32) \u\ <C\M{q,pC)u\ 

for all u G E_((7,p, ("). As p tends to zero, this implies that 

(D.33) \u\ <C\M{q,0)u\ 

for all u G E_((7,0, C) and |C| = 1 with 7 > 0. This means that M does 
satisfy the uniform Lopatinski condition for dx + iG'. □ 



E Appendix E. The viscous case 

In conclusion, we briefly discuss the case of a hyperbolic system with viscous 
regularization, that is, an x second order system with symbol 

Lu{p, r, = rid + A{p, + i^B{p, 

d d 

= rid + ^^jAj{p)s -ii^ ^ ijikBj^k{p)s 
(E.l) j=i j,k=i 

d 

j,k=i 

We make the assumptions of symmetrizability of the first-order part, i.e., 
existence of positive definite symmetric S such that SAj is symmetric for all 
j, and dissipativity of the second-order part, i.e. B{p,£^) = Ylj k^i^kSBj^^ > 
0, with no eigenvector of A{p,^) lying in the kernel of B(p,^) (the "genuine 
coupling" condition of |Kawj ) . 

The small- viscosity limit. In |MeZl[ inWMZ2l l(;WMZ31 K;MWZ4j . 

there was considered, under the additional assumption of constant multi- 
plicity of the first-order part, the problem of obtaining maximal estimates 



(E.2) 7hllL(M^) + Wm' < :^^ll/lli^(M+) + \9\' 

generalizing those of 1)4. 5() , where is the coefficient of viscosity ( "Reynolds 
number" ) in ()E.1|) and C = (r, ??) as in previous sections denotes Laplace- 
Fourier transform frequencies, for i/|d sufficiently small. Such small-frequency 
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estimates, together with analogous intermediate- and high-frequency esti- 
mates not depending on the constant-multiphcity assumption, were used to 
verify, respectively, existence of viscous boundary and shock layers in the 
ly ^ limit, with rigorous convergence to associated formal asymptotic 
series. 

A consequence of the matrix perturbation analysis of |MeZlj is that the 
reduced equations Cjjscous ^'^^ the viscous system for a hyperbolic mode 
associated with basis V satisfies the intuitively appealing relation 

(E.3) VGl,,,,^, = GV= {G,n,iscid + ^B{p, 0) V, 

where V = R{p, ^) is a matrix of right eigenvectors of A{p, ^) associated 
with the single eigenvalue r, and thus V*S = L{p,^) is a matrix of left 
eigenvectors, where S is the symmetrizer of the system. 

A well-known consequence of the dissipativity /genuine coupling assump- 
tion (see especially |Kaw| lKSh| ) is that LBR{p,^) = V*SBV is uniformly 
positive definite. From this fact, it is straightforward to see that, under 
the structural hypotheses of Assumption 16.151 the smooth hyperbolic sym- 
metrizer T, = V*S constructed in Section 6 near points of variable multiplic- 
ity serves also as a symmetrizer for the viscous system, yielding the desired 
maximal estimate (|E.2|) . This extends the results of the above papers to the 
variable-multiplicity case under A ssumvtion \6.1^ in particular yielding ex- 
istence of boundary and shock layers for MHD under the uniform Lopatinski 
(Evans function) condition. 

We remark that a continuity argument like that in Section 7 shows that 
the uniform Evans function condition reduces in the vanishing-magnetic field 
(H — > 0) limit to the corresponding condition on the limiting fluid-dynamical 
shock. The uniform Evans function condition is always satisfied for suffi- 
ciently small-amplitude shocks in either gas dynamics or MHD |PZ1 IFSj . 

Long-time stability of viscous shock waves. In |ZlUZ2llZ3llGWMZ2l 

there was considered under the constant-multiplicity assumption together 
with an additional technical hypothesis (satisfied for gas dynamics and 
MHD) that the glancing set have a certain "foliated structure" , the related 
problem of time-asymptotic n H'^ r\ stability, for s sufficiently 

large, of planar viscous shock profiles, for which the relevant symbol is ()E.H) 
with constant coefficients and v = \. Again, the constant-multiplicity hy- 
pothesis was used only in the small-frequency regime, this time to obtain 
maximal — > LP stability estimates, p > 2. Away from the glancing set, 
the only way in which the constant-multiplicity assumption was used was to 
show that the real part of "slow" , or "hyperbolic" eigenvalues was bounded 
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below by multiples of 7+|C| • But, it is a straightforward exercise to show that 
this is also implied by the existence of a symmetrizer SG^ > C^^{'~f + 
This extends the results of the above papers to the variable-multiplicity case 
under Assumption \6.15l in particular yielding long-time stability of Lax- 
and overcompressive type shock waves for MHD under the refined Lopatinski 
(Evans function) condition defined therein. 
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